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l. INTRODUCTION 


The abstract concept of orthogonality of functions (or other mathematical objects) 
Is a generalization of the notion of having two or more vectors perpendicular to one 
another. This concept arises naturally in a wide variety of physical and engineering 
disciplines. For example, the theory of orthogonal functions 1s central to the develop- 
ment of Fourier series and wavelets which are essential to signal processing. 

Classical Fourier series (real form) depend on the property that the trigonometric 
functions sine and cosine are orthogonal (on an appropriate real interval) in a formal 
sense that will be made precise later. As a consequence, a bounded periodic function 


f(x) of period 2r which satisfies the Dirichlet conditions! may be expressed in the form 


[+ a+) (a, cosax+b, sinnx) 


ne] 


where 


An 


x | (v cos n S 0 12 


b, = Pr pum npe 
are the classical Fourier coefficients. These formulas can be modified via a change of 
variable to accomodate any such function of period 2L. [Ref. 1: p. 529] 

This property can be used to generate other classes of orthogonal functions - 
polynomials, for example - that behave in very structured and useful ways such as in 
generalized Fourier series. In particular, specific families of these “classical orthogonal 
polynomials” have traditionally been used for solving problems arising in various areas 
of applied mathematics, physics, and enginecring, among others. 

This thesis develops the main ideas necessary for understanding the classical theory 
of orthogonal polynomials. Special emphasis is given to the Jacobi polynomials and to 


certain important subclasses and gencralizations. Much of the investigation will be 


I Dirchlet conditions: (1) In any period f(x) is continuous, except possibly for a finite number 
of jump discontinuities, (1) In any period f(x) has only a finite number of maxima and ininuna. 


made using the theory of hypergeometric power series and their q —extensions. The 
classes discussed in this thesis are but a small fraction of those identified and studied in 


the literature. 


A. CHEBYSHEV POLYNOMIALS 

The Chebyshev polynomials of the first kind, T(x), arise from an elementary trigono- 
metric consideration. As such, they satisfy various properties and identities which are 
easily derived directly from their definition, many of which are observable from their 
graphs (see below). This class of polynomials will serve as the model for some of the 
basic structure of more general classes. 


Foras 0 T2 ae chime 
T(x) =cos(n arccosx), -lsxsl 


Le. letting x cos Do Order. 

(1) T,( cos 0) = cos nô. 

Some immediate consequences of (1) are IT(x)| <1 for |x =< DPP Wl 
T,( cos <= ) =(—1), 0<k <n; in particular F(1)=1 and F(—-1)=(— I P ED 
Which can be seen graphically in Figure 1. 


l. Three-term Recurrence Relation/Differential Equation 


From (1), we have 
(2) Ty(x) 2-1 and Ti(x) — x, 
and by considering the identity 
(3) cos(a + b) + cos(a — b)= 2 cos a cos b 
with a = 10 b — 0, we obtain 
(4) Tal) = xT (x) — T,_,(x. 


Equation (4) is known as the three-term recurrence relation for T,(x) which to- 
gether with initial conditions (2) imply that 7,(x) is a polynomial of degree exactly x, 
called the n“ Chebyshev polynomial of the first kind. Note that the leading coeflicient 


of T(x) is 27! fortz 2. An inductive argument applied to this recursion siiommemenat 





Figure 1. Chebyshev Polynomials 


T,(x) is an even function if 7 is even, and odd if n is odd (see Figure 1). The first few are 


listed below: 
T(x) =1, T(x)=x, Ta(x)=2x7-1, T3(x) = 4x? —3x 
T(x) = 8x* — 8x7 +1, T(x) = 16x° — 20x? + Sx. 


Differentiating (1) twice with respect to @ yields the second order differential 


equation for T(x): 


(5) (1—x5)T,() — xT/(x) + n T, (x)= 0. 


2. Orthogonality of Chebyshev Polynomials 


Let 


I | ‘cos m0 cos nOdd. 
0 


If m z n, then using (3) with a 2 m0, b — n0 yields 


m+n mon me 


] | sin(m + n)0 sin(;/n — n)0 T 
ig as 


0 


(1 + cos 2a), 


j|- 


If m 2n 0, then by using the identity costa = 


T 


a m |? m PX | = SE 


0 
jí =, mO (on 
x 
loo -| dé = 7. 
0 
Hence, 
(6) | cos mÓ cos n0 dÜ — h;' oe 
0 
where 
pz) m=n=(0 
In = a2 oe 
and 


. Doo 
Onn = 


Ij 


is the Kronecker delta function. Changing variables via x = cos 0, we have 


dx —h,ló 


m,n* 


(7) | Trl) Tol) I 
E EU 


This important property 1s formalty known as the orthogonality relation for the 
Chebyshev polynomials. The reason for this terminology will become clear in the next 
chapter. 


3. Zeros of Chebyshev Polynomials 
| |.2k—] 


TT WJ 
2H ] Jr ; 


Setting  T,¡(cos0)=cosn0=0, we obtain 9=6,, 
=x, = cos0, ,, 1 < k < n. 

Thus all the zeros of T(x) are real, distinct, and may be regarded as the 
projections onto the interval ( —1,1) of the equally distributed points 0, , on the unit 
circle, as seen in Figure 2. Moreover, for 1 € k € n, an easy algebraic check verifies that 
DENMLSq v oaandsthereloke cx - < Xo s XeuruuQ: Hence, the zeros of T..(x) 
interlace with those of 7,(x). This interlacing of Zeros ts a striking feature of the plots 
in Figure 1. 

The zeros of Chebyshev polynomials, and of other orthogonal polynomials in 
general, are extremely important for applications to numerical analysis, electrostatics, 
and many other fields. 

4. Looking Ahead 

Many of the properties derived for the Chebyshev polynomials T (x) from their 
trigonometric definition (1), extend to more general classes of orthogonal polynomials 
via a general theory, elements of which will be developed in this thesis. Some of the 
many properties satisfied by these classes that we will derive include: 


l. Orthogonality with respect to a weight function 


t2 


Three-term recurrence relation 

3. Second order differential or difference equation 
4. Hypergeometric series expression 

5. Rodrigues’ formula 
6 


Generating function. 


The general approach we will take is to define these “classical orthogonal polynomials” 
via terminating hypergeometric power series, and from this prove (most of) the other 
properties. Howcver, because of this equivalence, manv authors choose to define a given 


class using one of these other characterizing properties. 


O - — - 


Figure 2. Zeros of Chebyshev Polynomials 


In order to understand the general theory, it is first necessary to define the ab- 
Stract concept of orthogonality in an appropriately defined “space” of functions. We 


turn our attention to these fundamental ideas in the next chapter. 


Il. BACKGROUND 


A. ELEMENTARY LINEAR ALGEBRA 
|l. Vector Spaces 
Let R° denote the collection of ali vectors (n-tuples), u 2 (a,, a,, ..., a,), where 
each ,eR,i=1,2,...,n. The standard Euclidean inner product (also referred to as the 


dot product ) of two such vectors u = (a,, a,, ..., a,) and v=(b,, b,,... , b,) 1s given by 


n 
<U y> = > ab, 
i=l 


The length, or norm, of a vector u e R’ is given by 


L 
2 
lul = /<u u> (Sita ji 


Two vectors u, v e R" are perpendicular, or orthogonal, if and only if <u, v> =0. 

The objective of this chapter is to extend these familiar notions to objects other 
than classical Euclidean vectors, in particular, the “vector space” of polynomials defined 
on a real interval [a, b]. 

A vector space V over a scalar field F (usually R or C ) ıs a nonempty set of 
objects called vectors, for which the operations of addition and scalar multiplication are 
defined. Addition is a rule for associating with cach pair of vectors u andy in V an 
element u + v, called the sum ofu andy. Scalar multiplication is a rule for associating 
with each scalar c in F and each vector u in V an element cu, called the scalar multiple 

ofu byc. [Ref. 2: p. 150) ` 

For all u,v, we V and c,de F, a vector space V must satisfy: 

|l. Additive closure. u ve V=u+ ve V 
Commutativity. u+v=v+u 
Associativity. u + (v + w) = (u + v) + w 


Additive identity. There exists a zero. vector, 0 & V, such that 0--u2 u--0— u. 


A 


Additive mverse. Kor each Wel, thete exists a vector —=ue Y, such that 
u+(=u=(—u+u=0. 


6. Multiplicative closure. ue V andce F = cu 8 V 


7. Distributivity. c(u + v) = cu + cv 
8. Distributivitv. (c+ du = cu + du 
9. Multiplicative associativity. c(du) = (cd)u 


10. Multiplicative identity. There exists a scalar ] e fF such that lu = u. 


Example I: R" (the model) 
Example 2: P,{a,b] = {polynomials of degree € N on the interval [a,5]) 
Example 3: P[a,^] 2 (polvnomuaals on the interval [a,5]] 
Example 4: C[a,b] = [continuous functions on the interval [a,b]) 
Note that Prla,b] < P[a,b] € Cla,b]. These vector spaces are sometime mics 
‘ferred to as function spaces. The interval [a,b] may be finite or infinite (em 
[a, co), (—oo, b], or (—oo, co) ) for our purposes. 

A subset U ofa vector space V is said to be a vector subspace of V if itis a 
Vector space in its own right. 

Example 5: P,[a,5] is a subspace of P[a,b], which in turn is a subspace of 
C[a,5]. 

Given a set of vectors (v, v, ..., v.) in a vector space V, atid Scala 
Ci, Cop ... , C, the vector c,V, + c,v, + - + cv, € V is said to be a linear combination of 
[V,, Y,, ... , VJ. The set of vectors [v,, v,, ..., V,) is said to be linearly dependent if there 


exist scalars ci, Ca ..., c, not all equal to zero, such that the linear combmanen 


ns 
CV, +¢,V, +--+, = 0. (Equivalently, at least one of the vectors v, can be expressed 
as a linear combination of the others.) Otherwise, the set (vj, v, ... , v,] 1s linearly inde- 
pendent. An infinite set S 2 (v,, v,,.., v, ... ] 1s defined to be linearly independent if every 
finite subset of S is linearly independent; otherwise S is linearly dependent [IXef. 3: p. $]. 
The vectors 4, v are said span Y if every vector ve V can be represented as 
a lincar combination of {Yi 6, ..., v. In this case, wesssrite I oV EUN "n v.}. The 
vectors (v,, v,, ... , v.) form a basis for V if they are linearly independent and span V. The 
dimension of V is the number of elements in anv basis. 

Example 6: The set {e,,e,,...,e,$ is the standard basis for R’, where 
e,=(0,0,..., 0, 1, 0, ... , 0) 1.e., the vector with a one in the 7 position amd Zerosvelse- 
where a 

Example 7: The set (1,x, x5, 5, x") is the standard bosistor z P p| IP m 
independence is ensured by the l'undamental Theorem of Algebra.) The dimension of 


P.[a,b} 1s therefore N+1. 


EIA eon A ee isethe standard basis for [a,b], and 
hence P{a,b] is an infinite-dimensional vector space. 
2. Inner Product Spaces 


An inner product on a real vector space V is a mapping 


<, >: V xV ->R 


such that for all u,v, we V and o, f e R, the following properties hold: 
l. Positive definiteness: <u,u> >0, and <u,u> =0 if and only if u= 0 
2. Svmmetrv: <u, v> = <v,u> 


3. Bilinearity: <au+fv,w> =a<u,w> + f<v,w> 


A vector space with an inner product is known as an inner product space. 
Example 9: V=R"; let constant "weights" w,— 0 be given, / 2 1,2, ..., n. 


Honu-iccscomangow-(b bsc b u, ve Y, 


n 
<U y> = U a; b; wj 
p 


If w,=1 for i=1,2,...,”, then this reduces to the standard Euclidean inner product, 
or dot product. Otherwise, this is referred to as a weighted inner product. 

The next two examples are commonly applied inner products on function space, 
and are analogues of the previous example. We assume a given weighs function 
w(x) > 0 in (a,b), integrable in the first case (e.g., continuous for [a,b] a finite interval). 


Pamplona Bla ol), Cleo 


b 
<f,8>=| Sogtomtoa 


Exampel: V= Plab] 


N 
<f g> = X Sogama) 
x=0 


(Positive definiteness is ensured by the Fundamental Theorem of Algebra.) 


The norm induced by the inner product 1s given by llull = <u, u> 2 


Example 12: For the inner products of Examples 10 and 11 therefore 


b 
ME (| | f(x) ie w(x) 4. ^ and 


N 
=D La JA 


x) 


respectively. These are sometimes referred to as ^L^-norms." 

Two vectors u, v e V are sald to be orthogonal, denoted u.Lv, if and only if 
<u,v> =0. The vectors u and v are said to be orthonormal | if ulv and 
lull = vi = 1. Note that the orthogonality of vectors in a space is determined by the in- 
ner product being used. 

The two examples which follow refer back to Chapter I, Section A.I. 

Example 13: Formula (6) shows that the functions ( 1, cos x, cos2x,... ) are 
orthogonal on [0, x] with respect to the uniform weight function w(x) 2 1. A similar 
computation shows that the same property holds on [—z, z] with respect to the weight 
Der 


function w(x) = + ! 


ee ==- ii PEN 


One advantage of preferring this inner product over the standard one lies in the com- 
putation of norms. Using |f] 2 / « f, f» , we have | H]| 2 2 and |] cos nx] =1 if 
nzl. Hence the functions ( 1/2, cosx, cos 2x,... ) are orthonormal on [—z, 7] with 
respect to the inner product above. Similar statements hold for the integral of a product 
of two sine functions on [—z, z], as well as for the product of a sine and a cosine. 
Example 14: By (7) the Chebyshev polvnomials {T,(x)} form an orthogonal 
class with respect to the inner product of Example 10 above on [ —1,1] with the weight 


function w(x) 2 (1 — x1)", 


2 Recall that «u,u» z 0. We remark that in the same way we defined inner product. earlier, 
it 15 possible to define a general norzm on a vector space which is not induced by an inner product. 


We remark here that Examples 10 and 11 can be unified into a single inner 


product on a “polynomial” space V via 


cee [se EE 


where da(x) is a positive Lebesgue-Stieltjes measure on a measurable set E possessing 
finite moments, Le, x" da(x) integrable, A — 0, 1, 2, .... In Example 10, E = [a,b] = R and 
da(x) — w(x)dx; the resulting expression is known as a continuous inner product, while 
in Example 11 the set E consists of a finite number of points (0, I, ... , N) S, and the 


associated measure gives rise to a discrete inner product. 


B. FOURIER SERIES 
Let ve V, and U be an n-dimensional subspace of V having some orthonormal basis 
{u,,...,u, }. (Any basis can be orthonormalized via the Gram-Schmidt process - see next 


section.) The vector v can be resolved into a sum of two vector components: 
(1) v=(v-w)+w 


muere we U and (v — w) E U. (See Figure 3.) The vector w is referred to as the 
orthogonal projection of v onto U. Since the vector (v — w) is orthogonal to every vector 


in U by construction, it follows that for each j= 1,2, ...,5, < v-w, u > =0, or 


(2) <y,y>=<w,u>. 


Moreover, since it lies in U, vector w can be expressed as some linear combination of 


a ul: 


Take the inner product of both sides with u, for cach j 7 1, 2, ..., ». From the assumption 


that «€ u,, 1, > = 0 unless ¿= J, we have the property that 


< W, u, > E UE u>. 


(3) c= < v, u; > 


11 


via (2) and the assumption that « uj, u, » =||u, | ^ I. 
Thus, 


n 


| 


and this vector represents the "best approximation" in U to ve V in the sense that of 

all vectors z e U, it is the projection vector w e U which uniquely minimizes the distance 

I v—z |. | 
Suppose now that U is an infinite-dimensional subspace of V (also infinite- 


dimensional), having orthonormal basis { u,,..., u,,... }. Then from (1) and (4) 


n 
v= (vw) +) < v, U, > Uu, 


i=1 


we may write 


(5) y = 


«y dcs Wy 


Ms 


ii 
— 


In the sense that 


(6) lim | v-w | = 0, 


le, the norm of the “residual vector” (and hence the vector itself) v-w — O as 
n > co. Formula (5) i$ known as the generalized Fourier series for v e V with respect 
to the orthonormal basis ( u, )Z,. The coefficients given in (3) are called the generalized 
Fourier coefficients of v e V. Statement (6) is Known as the norm convergence property 
of Fourier series, and the “minimization property” mentioned above extends to this 
infinite-dimensional case. 

Example 15: Let V=C[a,b], and {${x)}2, be an orthonormal basis of 
eigenfunctions (sonietimes referred to simply as an eigenbasis ) of V. Then f e V hasa 


Fourier series representation 


on 


fQ)- 2 ad) 


i=0 





U=spantu,,... 


Figure 3. Best Approximation 


with Fourier coefficients 


b 
(7 qe «f, d o m | f) ma) dx 


[n this function space context, norm convergence 


b 
lim | 
n-oo a 


n 


flx)— > abl) 


¿=0 








2 1/2 
dx = () 


Is referred to as mean square convergence, and is the least squares principle in regression 
analvsis. 
In particular, if V is equipped with the inner product of Example 13 and 


orthonormal basis 


Er COSO en een) 


on [—z,z] (see Example 13, Section A.2), then a suitable function fe V (and its 


2n-periodic extension on R ) has a classical Fourier series 


f(x) ~ + ag + > ( ap cos nx + b, sin nx) 


n=0 


where 


a,=<f, cosnx > => (anc OS ie arc 


b = <f, sinnx > => (x) Sin mx Zx. 


n 


I 
js 


as indicated in the Introduction. 


C. GRAM-SCHMIDT ORTHONORMALIZATION 

The Gram-Schmidt process orthonormalizes any set of linearly independent vectors 
in an inner product space. This method will be used in later sections for different inner 
products on the vector space P{a,b]. 

Begin with an inner product space V and any set of vectors (v,, V,, ... , V,, ... ), finite 
or infinite, such that any finite number of elements of this set are linearly independent. 


Recursively define a new set of vectors (u,, U,, ... , U,, ... } 


where ye =Y, We wit 


0 
3 By convention S a =O ANN E 


TEA si 
=] 


> 
i 
— 


w„=) <v,,u1?>U; 


w. 
— 


These new vectors (u,, u,, ..., u,, ... ) are orthonormal by construction and span the 
same space as the original vectors. Note that this process occurs in two stages: 
orthogonalization and normalization. The orthogonalization is accomplished bv sub- 
tracting w, the orthogonal projection of v, onto the subspace spanned bv 
{u,, U,,..., U,_,}. The component of v, which remains, denoted above as y,, is then 
orthogonal to the vectors {u,, u,,...,U,_,} aS shown in Figure 4. The normalization is 
then achieved by dividing y, by its norm, thus giving it unit “length”. 

|l. Legendre Polynomials 

Exanplew Ome been SA basis, lv. x, - } sand Uniform 
weight function w(x) = 1. The inner product is then given by < f, g > = J ! £@)g(x)dx. 
The Gram-Schmidt process yields the set 


X 


w px. = AW NET Ñ E = ( j 23 , an (s a 3! E 


as an orthonormal basis for P[—1,1]. Since this set is linearly independent, we can 
standardize the set by taking scalar multiples of these polynomials so that P(1) = 1. 


Members of the resulting orthogonal set 
oo dues | 3 
{P(x)};—-0 = 1 X F e — i. T [o — 3x), - \ 


are known as the Legendre polynomials on [—1,1]. If the normalized Legendre 
polynomials { p,(x) }@, are used as the orthogonal eigenbasis for a Fourier series, the re- 
sulting expansion is often referred to as a Legendre series representation; when 
Chebyshev polynomials are used, we obtain a Chebyshev series representation, etc. 

the Gran Schmidt process can always De used in this way to generate a class 
of orthogonal polynomials with respect to a given inner product (1.e., weight function) 
on a real interval, When using the Gram-Schmidt process from the basis 
SENS. n thevorthogonalization stage producing y, results in a set of monic 
polynomials, i.c., the leading coefficient of each polynomial is one. [n the normalization 
stage, we are dividing by the norm | y, || > 0. Thus the leading coefficient of polvnomuals 


in an orthogonal class is strictly-positive. In the next chapter, we will examine other 


span(u, ,u5) 





Figure 4. Orthonormalization 


ways to define these classes. It is the structure and applications of certain of these 


classes with which we will primarily be concerned. 


D. THE GAMMA FUNCTION 

The gamma function T(x) is a fundamental mathematical object that appears fre- 
quently in the representations of orthogonal polynomials as well as in many other ap- 
plications. This "special function" was developed as a generalizauon of the factorial 
function of the natural numbers. As we will see, the gamma function has the value 
(n — 1)! for the positive integers n but it is defined for noninteger values as well. 


A conventional definition for the gamma function is 


(8) I(x) = | eae tae O. 
0 


The positivity of x ensures that this improper integral converges. We now develop some 


fundamental propertics of the gamma function. Integration by parts in (8) yields 
(9) [ (x 4- 1) z xT (x). 


We now introduce the Pochhanuner symbol or shifted factorial, (a), , to simplify our 


notation. Forn>0, define 
(a), =a(at+ l)\(a+2)..(at+n—-—1), ifn>l 


rca), = 1. Letting a= 1 gives (D); = (026) ... (n) = n!. Note that for a negative inte- 
ger, (m), 2 0 if n2 m» O0. The shifted factorial can be defined for negative subscripts 
but we will not need this in our work with polvnomaals. 


Iteration of (9) n times yiclds 
(10) P(x +n) 5 G),EQO 


for every positive intcger n. Using this property, the gamma function can be extended 


to include negative rcal numbers by defining 


l 
(x), 





(11) B I(x+n) for —n«x«-—n- 1. 
Since this expression 1s undcfinecd when x is zero or a negative integer, the gamma 
function is not defined for those values. 

Letting x = 1 in (8) and computing directly, we have l'(1) 2 1. It then follows that 
[ (n - l) 2n! by letting x 2 1 in (10). Furthermore, 


oO = 2 
r( 3 ) - | UN gy = | e" du= Jn 
0 0 


Where the second integral can be evaluated by standard methods involving multiple in- 


tcgrals. 


Finally, we define a generalized binomial coefficient as follows. For x and æ non- 


x+a (x+a)! 
xl) au * 


negative integers, define 


For nonintegral a, define 





=) _ (a-- 1), a 


x (D, — T(x+1)T[(a+1) 


l. The Beta Function 
An integral related to the gamma function defines another useful function called 


the beta function which 1s given bv 


l 
(12) B(x, y) -| V = uw 
0 


for x , y » 0. We now establish an important connection between the beta and gamma 


functions. We start with an identity casily verified from (12): 
(13) B(x, y+ 1) = Bix, y) — B(x 1, y). 


Also from (12) 


1 


Bix + I, »-| (1 — Dn di 
0 


which when integrated by parts gives 
B(x + 1, y) = > B(x, De 


Substituting into (13), we obtain 





ar. 
Bi ya 


D(x, y) 2 —y 


which when iterated yields 


] 
Dp T NI E — | ee 


ur 


Changing variables from t to 4 , 


_ pn x-1/ı __L \n+y-1 
(14) B(x, ne | 1 Pa 


Taking the limit as n > oo, and using the fact that lim(1— c y =g" we have 


(15) B(x, y) = (x) lim e yh 
ZA nei 


(The fact that we can pass the limit through the integral on the right can be mathemat- 


ically justified.) If y 2 1, then (15) gives 


(x + 1), 


| X 
H. H 


BD: E(x) lim 


By direct evaluation using (12), 


l 
B(x,1) = | "PTS l | 
0 


Hence 
f (x + 1) 
LT) lim Zr, 
: noo nn 
which can be written as 
nn? 


(x) = lim 


noo x(x + 1), 


Noting that x(x + 1), = (x),., = (x).(x + n), we have 


x—1 





Lg nont 
POD = Te + n) 
which gives the form 
x-l 
16 (x) = lim —— 
(16) (x) = Jim EI 


(Equation (16) was Euler’s original definition of the gamma function. A separate "esti- 
mation” argument may be used to show that this hmit mathematically exists.) Thus by 
(15), 


(v T 
! X+y—1 
B(x, y) =T (x) lim AAA 
(x, y) = T (x) lim e 
n! n? 71 
Then by (16), we have the useful identity 
[ (x)T (y) 
(17) B(x, y) m [ (x + y) : 


We will find this identity useful in understanding the Jacobi polynomials in 


Chapter IV, where it becomes necessary to evaluate a related integral: 


| | UI Mi dx. 
-i 


We remark here for future reference that the formal change of variable x — 1—2t can be 


used to transform this integral into 


1 
a | "(-N d= 2 Blatt, B41) 
0 


repair Clot) P+) 


x [(@+ß+2) 


ll. GENERAL THEORY OF CLASSICAL ORTHOGONAL 
POLYNOMIALS 


In this chapter, we examine some of the characteristic properties associated with 
classes of orthogonal polynomials. Some of these properties provide alternate means 
of defining a class. These alternate definitions often provide a straightforward way of 
producing a specific result that may be very difficult to derive otherwise. 

Throughout this chapter, we let { p,(-x)}, denote a set of real polynomials with 
bo lol degree n, 1.e., 

n n—1 
(1) D ac Ss N. 
Recall that these polynomials are said to be orthogonal on an interval [a,b] with respect 


to a continuous weight function w(x) > 0 on (a,b) i£ 


b 
(2) < Pm Pa > = | PREMIER = hy" Syn 
a 


where the normalization h-' 4 0 is chosen to simplify the expression of certain formulas. 


iNote that since 


b 
Lo | Los Gown m n. 


a 


it follows that h, > 0. 


A. POLYNOMIAL EXPANSIONS 
We begin by showing that any real polynomial q,(x) of degree m on [a,b] can be 


written as a linear combination of orthogonal polynomials { p,(.x)}"_, : 


(3) aml) = Dim PAX) 


m 
i=0 


for constants « 


“tym 9 


The proof is by induction on the degree mt. Since g (x) is a real polynomial. we 


A ON 


write 


21 


A AS M E M 


where a, #0. For m=0, (3) reduces to @ = a9 X using the form for p(x) given in (1) 
and &,, Is uniquely determined. 
For the induction hypothesis, suppose that for m > 1, we can write any polynomial 


of degree m — 1 as a linear combination of Í p,(x) 155 : 


m—1 
Qm 1c) = > Xi m-i p(x). 
i=0 
Since q,(x) — (a,]k,)p,(x) 2 q,.(x) is a polynomial of degree m — 1 the induction hy- 


pothesis implies there is a representation 


m-1 


a 
Inl) (| pal) = >, it PA) 
m i=0 
Now set &„m = (a„/k„) and the result (3) follows. [Ref. 4: p. 33] 
Using the theory of Fourier series developed earlier, we next determine the coeffi- 
cients a, explicitly. For £ — 0, 17 mete sn hanoi 


on» 
Oe)” oie T 


Then by construction, { 6x) }%, is an orthonormal set of polynomials. Writing (3) as 


Iml*) =) Sim PAX), 


m 
i=0 


we see that the right-hand side can be interpreted as a (terminating) Fourier expansion 
of q,(x). Hence the results of Example [5 in Chapter II, Section B may be applied. In 


particular, by (7) in that section, the coefficients c,,, are given by 


b 
Cim=< In: O> = | qux) ix) w(x) dx. 


d 


Changing back to the old variables, 


b 
aim = h; oe S ^] q^ x) pix) w(x) dx. 


a 


We are now in a position to show that the orthogonality property 


b 
(4) | p. (x)p,(x)w(x)dx =Ü ,m=n 
a 


can be expressed equivalentlv as, 


b 
(5) | P.) x^ w(x)dx 2 0, m «n. 
a 


To see (5), substitute the form (1) for p,(x) into (4) where m « n. The lineanty of the 
integral gives (5). On the other hand, since x" is a polynomial, we can write x” as a 
linear combination of the orthogonal polynomials, so (5) gives (4). Note that (5) implies 


each p,(x) is orthogonal to every polynomial of lower degree. [Ref. 4: pp. 33-34] 


B. THREE-TERM RECURRENCE RELATION 
The three-term recurrence relation is a useful result which holds for any three con- 


secutive orthogonal polynomials: 
(6) py) — (As x *- Bp, AQ) — Copa) n 2,3,4,... 
where 4, , B, , and C, are constants given by 


Kn Sa Sn-1 
A, = k > 0 ; B = Al k, = k, : C, = 


ky, k 2 An? 





Wiesrecttrence relation is valid for n= 1 if p_,=0 with C, arbitrary. In this case, the 
formula for A, also holds for n= 1. (In the contrapositive form, this statement 1s a 
powerful tool for showing that a polynomial set is not orthogonal.) [Ref. 5: p. 23- 

To prove this, we begin by considering p,(x) — (&,/K, ,) x p, (x), a polynomial of de- 
gree no greater than (n — Il). We expand it in terms of the orthogonal polynomials 


(p(x))z, via the technique of the previous section to obtain 


23 


n—| 


ky 
(7) Pals) = FA x pat) m 2 tui p. 


i=0 





The coefficients oe, arerdetermined by 


in 





k 
(3) "ES TT h, k 3 a , pix) > muore 
n=l 


Because 
b 
€x p, (x). p(x)? — | x p, (x) p(x) w(x) dx 2 € p, (x) x p(x)? , 
it follows that for i < n — 3, x p(x) is of degree no greater than (n — 2). Hence 
<p,-i(x) , x p(x)> = <xp,-i(x) , p (x)> =0 ,1<n— 3, 


since p, (x) 1s orthogonal to every polynomial of lesser degree. Thus the constants 


Ario a nn. are all zero leap oc ETE: . With this knowledge, (7) 





a 
becomes 
ky 
Ppa ox TE XPa) S aa a n Pn a) a S p n 
€ 
Setting 4, 2 kk, ,, B,20, 4,4, andC,2—28, ,, , then rearranging terms gives (6). 


To determine C, explicitlv, we write 
C, = ze An-2,n-1 m A e , Pn-2(*) a 


from (8). Since 


n=] 


xpo E 


— n—2 zd ... 
p K [kni B ] 


n—2 
—- m T DUUM no 


A 





we can Write 


= 
A 
C, m harm Spb Pre) +) sca rd 
n= —=0 





Using the properties of the inner product and the orthogonality property of the 


polynomials, we obtain 


A, h 
C, = h,—2 pu SERO) Spo = m c» s CIEN 


astelven in (6). [Rel 6: p. 8] 

Into equation (6) we substitute the expanded forms of the polynomials 
DIES»). px) from (1) and the constants A,, C, from above. [quating coeffi- 
cients of x' gives B,. Since k,>Oandh,>0, it follows that 4A, > 0 and C, > 0 and the 
proof is complete. 

The nonnegativity of the constants A, and C, is important for the converse of the 
result in (6). Favard showed that the existence of a three-term recurrence relation in the 
form of (6) implies that the polynomials of the set are orthogonal with respect to some 
weight function over some interval using Stieltjes integration [Ref. 7]. 

We observe that in order to generate the polynomials of an orthogonal class with 
the recurrence relation (6), we need the sequences of constants 4,, B,, and C, together 
with two of three consecutive polynomuals in the class. Other techniques provide what 
some authors call a pure recurrence relation. requiring only two of three consecutive 
polynomials to define the class, because the constants as functions of n are contained 
explicitly in the recurrence relation. The recurrence relation derived in Chapter I, Sec- 


tion A.l for the Chebyshev polynomials is an example. 


C. CHRISTOFFEL-DARBOUX FORMULA 
The Christoffel-Darboux formula is an important identity which can be derived from 


the three-term recurrence relation. The identity 1s 


Bor) 


- h, k, 
(9) han) = 2 = 


j=0 n+] 
To prove (9), note that from (6) we have 


Pr) (Art Bde) — Cpr) 
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which we rearrange to give 


Bis C; 


l j j*1 
(10) xp(x)—— Dux) EE EE 
pi A J+1 Ay. @ DLE l 


and similarly 
10 D; +1 C 
(11) YP S m (D — AO) Puer S 
Aui Aj J+1 
These recurrence relations are valid for j 2 0 if we set C, 2 0 and p (x) z pay 


Multiply (10) through by p(y) and (11) by p(x), subtract the results, and then mul- 
tiply through by #, to obtain 


h; k; 
h =A AA A= TE LAA — aA] 
J 


hy ky 
ur IS [ pj COp rata to). 


Summing over j from 0 to n yields a telescoping series 


h. k 
(sy) 2» iy P(x)p{ y) = e PCO TE ay 


from which the identity (9) follows. [Ref. 4: p. 39] 
Now subtract and add the quantity p,(x) p,.,(x) to the numerator of the nght-hand 
side of (9) and let y tend to x to obtain a limiting case of the Christoffel-Darboux for- 


mula: 


3 





Kn 
| Pa (2) Pac) =p ea Je 
j=0 n+] 


We will use (12) in the next section. 


D. ZEROS OF ORTHOGONAL POLYNOMIALS 
In Chapter I, we observed that the zeros of the Chebyshev polynomials {7,(.x)} are 
real, distinct, and lie in the interval ( —1,1). This principle extends to any class of 


polynomials orthogonal on the real interval [a,b]. 
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To see this, choose 5» 2» 0 and suppose that p,(x) is of constant sign in [a,b]. Then 
< p(x), p(x) > #9, which contradicts the assumed orthogonality. Thus by continuity 
(and the Intermediate Value Theorem) there exists a zero x, € (a,b). 

Suppose that x, is a double root. Then p,(x)/(x — x,? would be a polynomial of de- 


gree (n — 2) and so 
2 j oe) > > 0 


which is a contradiction. Thus the zeros are simple. 


Now suppose that p,(x) has exactly j zeros x,, x, ... , x, e (a,b). Then 
E e (u= x;) E ej (Eo ee E > mr E 

where q, , (x) does not change sign in (a,b) and 

SH) ik — xx — x) — xj)? mdi (x), (x — xi (x = x)” ... (x — >) = 
Since 

r (S) nnm a 
and 
p. ee O ne N ee 


then it must be that j > n. The Fundamental Theorem of Algebra precludes j > n and so 
we conclude j 2 n. [Ref. 5: p. 236] 
Thus all the zeros of p,(x) are real, simple, and lie in the interval (a,b), and so may 


be ordered 


yr i c cce D. 


The interlacing of zeros of p,(x) and p,.,(x) follows from the Christoffel- Darboux for- 


mula. Recalling that the leading coefficient k, is positive for all p,(x), then (12) gives 
(13) e IRA Cp 0 aaa. 
Bern and y be adjacent zeros of p.(x). Then 


(14) (O 
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since the zeros are simple. At these zeros, inequality (13) reduces to 
—p'n(u) py au) 7 0 

and 
pP’) Pr) >. 

Multiply these two inequalities together with (14) to conclude 
Pha 0 Pu) « 0, 


so p,,,(x) has a zero between each pair of consecutive zeros of p.(x). 
Now let x,, denote the largest zero of p,(x). Observing that p,(x) — oo as x — 00, 


we must have p',(x,,) > 0, and so by (13) 
DESC) <0 


But p. (x) > co as x > oo, SO p,,,(x) must have a Zero to the right of x,,. Similarly, 
P.. (X) must have a zero to the left of xin the smallest zero of p(x). Thus all n+1 zeros 


of p,,,(x) are accounted for and interlace those of p,(x) [Ref. 8]: 


a < X] x+1 < ln < X2 nl < ae S Xp nel = Xen < Xk+1,n+1 < Pu < Xan+l < Xan < Xn4-1,n7-1 < b. 


E. GENERATING FUNCTIONS 


The function f(x,t) having a formal power series expansion in t 


Fix) = ) f 9r" 
n=0 


— 


is said to be a generating function for the set { f,(x)} [Ref. 9: p. 129]. For an appropriately 
chosen generating function, the generated set | f,(x)} is a class of orthogonal polynomials. 
By defining a class in this way, properties of the polynomials can be derived from the 
generating function itself. 


For example, consider 


(15) Fix) = (1 - Zer + PY =) flo" 


n=0 
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which for fixed x is the Taylor series of P(x,t) centered at t = 0. I( we restrict x to the 
interval [ —1,1], then by considering the singular points of F(x,t), we may conclude that 
Hoeeseties 1s convergent for |r| <1 (Ref. 4: p. 25]. Thus we can determine /(x) for 


D= 1.2... oy 





1 a = 
(16) eS P Mox Ete 
Equation (16) yields f(x) 2 1, f(x) x, f(x) ase, etc. We also note from 
(16) that 
lence = 
fl) =F pr A Y Jee 
l EE 
Ser (1-10) a 
= |. 


We now derive some basic properties of the f(x) from the generating function in (15). 
l. Recurrence Relation 


Differentiating (15), we find 





Kr E NR N pen 
(17) = (1 — 2xr +) 22 = Ah (xt 
CE 2:23 ` n—1 
(19) < =(x— D(1 — 2xr + Ë) = 2 nto | 
Simce (x — 1) oF - mu — 0, we have 
Ox Ot 
(x — )? f, Cor" — 2 n f(x)! = 
n=} n=l 
which becomes 
» x f(x)” — n f (x) = DO 
nzsl n=1 n=] 
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Since f(x) = 1, we can start the sum on the right side at n — 0, then re-index so that it 


starts at n = ] again. We then find 


S [oce cm SR MN J-Y ao 
=! 


Zl 


— 


Equating coefficients of 7" then gives 
(19) x fn (x) - nfe(x) - f(x), nz lH 


Rewrite (17) as 


(20) (17 2x 2)? 2. pog , 10. 
n=l 
Substituting the appropriate expressions from (20), (18), and (15), respectively, into the 
identity 
(1 — 2)(1 — 2xt + £0) ?? — (20 (x - 9(1  2xt NE 2 (1— 2t (0). 
now gives 
(1-5)? fro"! — Qr D a - Ys I eor 

nzl n=} 

which, upon rearranging, becomes 
> fs (9 En | Q9) yt Y anf) y ds xr. 
n=l 
Equating coefficients of t" and gathering terms yields 
(21) (2n + Ds) m fa GO) fa O), n>. 
Substituting (19) into (21) gives 


x Pa = fai (x) = (n a DÍ 


which by a shift of index from n > n-— 1 becomes 
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(22) xfi (0 m Om nf a). n>2. 
Substituting again from (19) gives 
(23) (x^ — 1) f(x) = nx file) — n £,_ (x). 
Multiplying (19) through by (x? — 1), we have 
x(x? — 1) fi'(x) a DA = (<° — 1) £, (O. 
Now substitute for (x? — 1) f’(x) and (x? — 1) f,-ı (x) from (23) and gather terms to get 
(24) nf o) - Qn- D xf a) - (r— Df), n2 2. 


Equation (24) is a tAree-term recurrence relation for { f(x)}. The advantage of 
this form is that beginning with f(x) and f(x), we can now generate any member of 
{ f.(x)} by iterating (24) and thus avoid the differentiation in (16). Note that since 
fx) 91, f(x) 9 x, (24) implies that [ £.(x)j is a set of polynomials. [Ref. 9: pp.159-160] 

2. Ordinary Differential Equation 

We continue the same line of reasoning to extract additional information about 

{ f,(x)}. Differentiating (22) yields 


(25) ia == (+ Df (o) 


IcmEUI9) we have f(x) and, after diflerentiating, £, (). Substituting these ex- 


pressions into (25) gives 

x [x fs (x) — Qr 1) Q9)] 9 Q0) — GC DEx A G0) — ne] 
which when rearranged becomes a second order ordinary differential equation 
(26) (1 _ FACES — 2x f. (x) + n(n + 1) f(x) = 0. 


The [ f(x) are solutions of (26) for n — 0, 1, 2, ... . (Ref. 9: pp.160-161] 
3. Orthogonality 


Rewriting (26) as 


(27) —— [(1— x) e | + a + 0/0) 50 


al 


we recognize the structure of a Sturm-Liouville eigenvalue problem. Since the points 
x= + 1 are singular points, we require that fí(x) and f. '(x) be finite as x > + 1. From the 
associated theory of the singular Sturm-Liouville problem, we conclude that the { f/(x)] 
are orthogonal on the interval [ —1,1] with weight function w(x) = 1. 


Alternately, we combine (26) with 
d 
LE (1-2) fy'(x)] + mln 045,6) = 0 


where n + m, to obtain 


falo E [1-4 00] 2 500 A [Q1 9509] 
T [n(n 4 1) — m(m t 1)] fA(x) f(x) 2 0. 


(28) 
Since 


- [(1 - 3,560409 —/562/59] 
- [cof - 37) 5G) OO 0) | 
LIE] e [0 =") nl} 


we can write (28) as 


(nO Sat (2) Lo 00401] +10? m4 nm] fate) ft) =0 


Or 
(n = m(n + m+ DAOS = [O — x) 509/509 — 45094509) . 
Integrating from x = —1 to x = 1, we obtain 


l 
(n — m)(n + m + 1) | Srl) fet o [(1 — x?)15/ 09,569 750945691 |: 
_! 


Ssincel(l = 2) = 0 at y = + we ha 


l 
(n — m)(n + m + j| DNE v =m 
-1 


Recalling that n + m, we conclude 


| e) ee =o 
=! 


i.e., the polynomials {f(x)} are orthogonal on the interval [ —1,1] with respect to the 
weight function w(x) = 1. [Ref. 9: pp. 173-174] 

These are the Legendre polynomials that were previously defined using the 
Gram-Schmidt process. Thus the f(x) of (15) are in fact P,(x) and (15) may be written 


Oo 


(29) (1 — 2xt + mi = > neo 


n=0 


establishing the equivalence of the generating function definition with the Gram-Schmidt 
definition. 
Legendre and Laplace concluded that the P(x) in (29) were polynomials of de- 


green in the variable x by examining a series expansion of the function 


oo [n/2] 


(1-2xr r^) 7 E » r” > CD Om UT nm (omen 
0 


mom yl 
ae. ie im! (n— 2n)! 


They reasoned the orthogonality directly. From (29) we write 
(1-2xr+4r 7 = > P(x) r 
i=0 
and 
(1—2xs+s2) 12 = » P (x) 5. 
j=0 
Multiplying these power series together via the Cauchy product formula vields 
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asus T9 k n-k 
T I-2xr-r? /1-2xs+s? hr NEN! 
| 2d ht | id. 


n=0 L i+j=n 


Integrating from -1 to 1 with respect to x gives 


l E l 
| in | Pix) Pix) is |e 
/ 2 2 4 " 
= ] -2xr cr Ji — 2xs s n=0 i+j=n l 
Through tedious calculation, the left-hand side of this expression becomes 


| L + rs 


nr 


E I= Jrs ’ 





a function of the product rs. From this we conclude that the coefficients of the terms 


in the series on the right-hand side are zero when / and j differ, 1.e., 


l 
| P(x) P(x) dx 20, ij 
_1 


and the orthogonality is established. Finally, 





(1-2xr- 1)? = = 





ON 
n 
=) 
n= 


or P.(1) 2 l and so the P,(x) are in fact Legendre polynomials. [Ref. 8] 
In a similar fashion, the norm of the Legendre polynomials can also be obtained 


from the generating function. Let 


l 
s mA m | [neo a 
=i 
Square both sides of (29) and integrate with respect to c on [ —1,1] : 
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Y y J. COPI CX) ás) ass -| ee ke 


2 
m=0 n=0 a —] l = 2 di +! 


M 


By orthogonality, the left-hand side is zero unless mm =n, while the right-hand side is 


integrable in closed form: 


l 





2.6 "= ET In (1 — 2xt ^ (^) 


IU gens 
ol jp 


This function can be expressed as a difference of two logarithms, each of which has a 


—1 


convergent Maclaurin series expansion in ( —1,1). When combined this yields 


2,5 =) TP je 


n=0 


Comparing coefficients of 1 on both sides gives c, = Sr 

Since there is no systematic theory for determining generating functions, finding 
one for a polynomial class can be a problem. The work above bears this out. Unfor- 
tunatcly, the proofs above do not easily generalize to related classes. With this in mind, 
let us summarize the key steps in proving that (15) is a generating function for the 
Legendre polynomials. First we cstablished that f(1) =1 for 1 >0. Next, we showed 


Mat 


for f(x) a polynomial of degree n in the variable x. Finally, we showed that these 


polynomials were orthogonal on the interval [ —1,1] 


[ DU COE Den 
Si 


These three points are sufficient to show that the generated class is the Legendre class 


of orthogonal polynomials. 
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We now provide an alternate proof attributed to IIermite. Our interest in this 
proof is mainly the technique which suggests a method of generalization that we will 
take advantage of in the next chapter. 


We begin with (15) reproduced below 


(1-28 + PR) ar. 
n=0 
Multiply both sides of this equation by x^and integrate from -1 to | to obtain 


x“ dx 


(30) [ —— E 2 [ x* f (x) dx. 
-ı  /1-2x+t mo oH 
Now change variables from x toy via 
(1—-2xr i) "^1 ty 
giving 
al 


uc c = - y, dx —(l — ty)dy. 


The left-hand side of (30) becomes 


f E (1-7) + | 


(1 — 1y) ay 


or simply 


J E r (1 — y) +| 2 


Expanding the integrand of this last expression via the Binomial Theorem, we obtain 


k 
NT 
AAN 
Os 
JS a 
= 0 
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When written as 


Y E p i (1-3 y^ 4) T, 


J= 


we can identify the form as a polvnomial of degree k in the variable t. Comparing this 


form to the right-hand side of (30), we conclude that 


| x* f(x) dx 20 , forn>k, 


1. €., f(x) is orthogonal to all polynomials of lower degree. The same argument as before 


Bes (1) — 1 for 7 >0 and the proof is complete. 


F. HYPERGEOMETRIC SERIES 
The term “hypergeometric” was used in 1655 to distinguish a series that was “be- 
yond’ the ordinary geometric series ] 4- x -- x? - c. In 1812, Gauss presented the power 

scries 
ab x  ala+ Db(b+1) «e? — a(a- D(a-- 2)*(b -- D(b --2) y? 


B 07 NOE oe cle + Ile + 2) 5 





c x 0, —1, —2, ... which is known as Gauss' serics or the ordinary hypergeometric series 
[Ref. 10]. 

Convergence of this serics for |x| « 1 follows directly from the Ratio test. By 
Raabe's test, convergence can be shown for [xl 2 1 when (c — a — b) 2 O[Ref. 11: p. 5]. 
Gauss also introduced the notation j,Z[a5;c;x] for this series. Note that 
iab ;c;x] may be considered as much a function of four variables as a series in x. 
MS! 12:p.1] 


With the shifted factorial, the ordinary hypergeometric scrics can be expressed 


owe 


n=0 


E 
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Delow are some examples of important functions which can be expressed as ordinary 
hypergeometric series. 

Example] eeltm 9 Ze me xl 

Example 2 sn Q Sx A I 

Example 3a) =: AR SA 

The generalized hypergeometric series ıs formed by extending the number of pa- 


rameters, an idea attributed to Clausen [Ref. 11: p. 40]. 


(a,)4(42), = (a), n 


ELA a a Dl TURNER EMEN (PS: i 


nz 


Note that since (a),.,/(a), 2 n + a, a hypergeometric series = c, x^ is characterized bye 
fact that the ratio c,,,/c, of coefficients is a quotient of two polynomials in the indem 
1.e., a rational function of n. 

The Ratio test can be used to show convergence for all values of x when r € s and 
for |x| < 1 when r=s+ 1. When ros 1, the series diverges for all x #0 angie 
function is defined only if the series terminates. The series terminates when one or more 
of the numerator parameters a, is zero or a negative integer [Ref. 11: p. 45]. This is an 
important characteristic of the hvpergeometric series that will be used later. A power 
series that terminates gives a polynomial which is defined for all x. In this case, the pa- 
rameters b,,..., 6, may be negative integers as long as the series terminates before a zero 
is introduced into a denominator term. 

Examples of the generalized hypergeometric series include familiar functions such 
as: l 

Example 4: (1 + x 2,F[ —a; — ; —x] 

Example 5: “= ,F,[ — ; — ;x] 

Example 6: sinx =x ,f,[ — ;3/2; — 2/4] 

Example 7: cosx = F| al 2 eee on 


Example 8: The Bessel function of order « 


(x/2)°9Fi[ — 3a +13 —x7/4] 


7,18) = r(x + 1) 


X 


where dashes indicate the absence of parameters, ie-, when r=00rs=0. We ws 


use a common alternate notation : 
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CRT 
5 l nr x | 
L D. 
for either the ordinary or generalized hypergeometric series. [Ref. 12: p. 4] 


l. Chu-Vandermonde Sum 


The Chu-Vandermonde sum 
4 L] i s n 
Fi —n,a;c; 1] = 
is one of many useful summation formulas. Since this one will be used in a later chapter, 


the proof is provided below. 


Basically, this is a consequence of the General Binomial theorem 


(y 
k=0 


Starting with the identity (1 — x)-*(1 — x) = (1 — x)-*-*, expand both sides. Us- 
ing the Cauchy product on the left side and the Gencral Binomial theorem on the nght, 


we have 


s.) ETÀ. — E 


nz n 


eS 


—k . . 
. Equating coefficients of x^, we have 


n 


>. (a) (P), y E (O) 


(Ada (1), 


k=0 


In order to express the left side as an ordinary hypergeomictric series, ,/,, mul- 
tiply both sides by (1), and use the identity (1),/(1),_, = (—I)*( —”), to obtain 


31) B OS iar Gs bs 


c (l); 


NEXT, 


3? 


( (5), ., = (DUDO k 
=-(- - DDr mm 


[ (—b—n+ 1), 
, (—b — n + l), 
gie rere 


Using the above result, equation (31) becomes 

` mika L R n 

E (Me (bnth (-b-n-D,- 
Let c= —b—n+ 1, and substitute to get 


=c+a=n+1), 
Ema ee 


(c), 
u (-1)"(—c+a-—n#l)...(-—c+a) 
y (c), 
B (c— ac n— l0)...(c — a) 
7 (c), 
>= een 
POCHE 
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IV. JACOBI POLYNOMIALS AND SPECIAL CASES 


This chapter focuses on the orthogonal class known as Jacobi polynomials, P(x). 
Until recently, the classes of orthogonal polynomials considered "classical were usually 
given to be Jacobi, Gegenbauer (also called ultraspherical), Chebyshev (of first and sec- 
ond kind), Legendre (also called spherical), Laguerre, and Hermite. The Jacobi 
polynomials hold a key position in this list since the remaining classes can be viewed as 
special or limiting cases of this class. Today, the classical orthogonal polynomials are 
taken to be special or limiting cases of either of two very general orthogonal classes 
known as the Askey-Wilson polynomials and the q- Racah polynomials, between which 
we will establish a formal equivalence. Because of their complexity, these classes are 
described in Chapter VI after the necessary additional theory has been developed. 

The results derived in the text that follows are arranged in tabular form by class at 


the end of this chapter. 


A. JACOBI POLYNOMIALS 
l. Definition / Orthogonality 
The Jacobi polynomials P" x) are generated by applying the orthogonalization 
step of the Gram-Schmidt process to the standard basis (1, x, x?, ... ) of P[ —1,1], with 


respect to the weight function given by a continuous beta distribution on [ —1,1] 
w(x;a, f) e (1 RUHR 


for & > —1, 8 > —1, i.e., 
l 
| Bo alana) s O O ou 
_1 


Ihe Jacobi polynomials can also be represented by hypergeometric series 


(1) pe Bye = e+ Da F, [^ H + w + f + | 


$ 
n! a+] 2 


This set of polynomials is thus standardized (as was done for the Legendre class in 


Chapter II, Section B): 
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H n! 


us = (oer D I 


We shall demonstrate that these two characterizations of the Jacobi polynomials are 
indeed equivalent. 

To show that the polynomials defined in (1) are in fact orthogonal with respect 
to w(x;a, f) 2 (1— x)'(1-- xf on [ —1,1], it suffices to show that Pf") is orthogonal 
to one polynomial of degree m for m=0,1,2,..,n— 1. This is because anv PÜ”(x), 
0 X m € n — l1,can be expressed as a linear combination of such polynomials. While any 
polynomial of degree m could be used, we choose (1 + x)” for reasons that will become 
apparent. 


To establish orthogonality, we consider 


l 
< pe B) (1 +x)" > -| pe Po (1 20" (1-0 (1+ 0 dx. 
—| 


By (1), the right-hand side becomes 


n 


(2 + l), y A | (Ae tx 


m (a + 1); k! 2“ =i 


Using the last result of Chapter II, Section C.1, the change of variable x = 1 — 2t vields 


(a+ 1), i Í ( —n), (n + w + f + 1), 


m Pe" m 1 DOSES NIIS EE D 
° a ki 


k=0 


E (2a 1), y ( sima + F(k 4 & - D Fn 4 f + HE) 
on! k! ) («+l),Uk+m+a+ß+2) | 


k=0 \ 


Identities for the gamma function allow us to simplify this expression to the form 
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n 
TA Ia +n+ 1)T( + f + 1) ( —n), (n + x + f) + 1), 


[ (m 4 a 4 f +2) (m+ a+ Ü + 2), k! 
KA 


which can then be written 


ee eg ee Ge Beal) [na BRL 
[D (m 4 w + f + 2) a m+ x+ f + 2 


The Chu-Vandermonde sum from Chapter II allows us to write this as 


ett mt hi) _ (tian, 


D (m 4- a 4- ff 4-2) (m+a+ß+2, 
Thus, 
C,, mn 
pr one, > -h. m=0,1,2,...,n—1 
where 


A A ae BD 


n= Bp cS 


which justifies the orthogonality. It is possible to extract the value of [A^ P] by modi- 
fying this argument, but we defer this computation till the next section, where it will be 
easier. 
2. Ordinary Differential Equation / Rodrigues’ Formula / Norm 
We begin deriving the ordinary differential equation for the Jacobi polynomials 


by noting the general formal result: 
emery TET _[at+tl,..,aq4+1 
(2) SP | cm d 2d l r E 
dx Dare O Di cess b + 1,..., b. + | 


This can be seen from the definition of the generalized hypergeometric series by 





writing 
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o ED (ak... (a), xt 
de G A dx M4 (bi)y (be K 
p Die (ay xem 
(b). (b)y (k= 1)! 
(ar (a) xt 
k=0 (b, Jeti (Oi KG 


Noting that (a),,, 2 a(a t 1),, we have 


oo 


7 B "n | LE T (a D,..(a + Do x* 
ERE zh x = Oa e sa m 
a b. bb, £a (D, OD, A 





from which (2) follows. Differentiation of the hypergeometric series is justified by re- 
calling that the power series is convergent for | x] <1 whenr=s+l. 


We apply this result to the Jacobi polynomials to obtain 


+ | -un+ta+ß+l ı_ 
Lo EA j P O p A |) 


dx (à a+] 
TIN A a a : -n+ l,n+a + p t NE 
E n! a+ | ES c EMO ' Ee 
_ (m+ EB EIC mmn =; peti. BHD) 
2 (n — 1)! (a + ia, n—l 


which simplifies to 


TREE CG E 
(3) PA) = TATE pet PD 


With these results established, we now consider the orthogonality property of 


the class; eser 
l 

(4) IP | pe Pc) m exem [AP PT? 5, 
E i 


By (3), we can write 
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l 
C NER RES (x, ff) y B| d p-i, f—1) 
Inn = "ud | P esce) E ei E J ax: 


Integration by parts yields 


o E EO- ) | pest ads 


noting that the boundary terms vanish, so that 


(6) d = m Po ps BM) q,+1 (x) (! — x)” (1 4 x7 ax 
where 
( — eu) - (1-x)) a P&P) — «(14 x) PP) 4 g (1 — x) Pm P(x) 


is a polynomial of degree (n + 1). Thus by (6) we may express g,.,(x) as a linear combi- 


nation of Jacobi polynomials 


n+! 
(2—1, —1) 
Ing 0) = p c; Pj Š (x). 
j=0 
We would now like to show that c, = c, = ::: = c, = 0, so that only the last term 


of the summation survives. The constants c,, j — 0, 1,2, ...,n can be determined by 


substituting this expression for q,.,(x) into (6) and using (4). 


je A) x -| p PX (x) po P ) (1 — x) (1 + x) dx 
n+l 


2? q n n n s E GE 
= pope i | p (x) P, (x) (1 = x) (1 e x) Ue 


For each m € n — 1, by (4) the left-hand side 7, , = 0, but the right-hand side is 


zero unless j= m + l <n. Thus 
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= (2-1, B=1)J]-1 + 
(zZ Cm+1 Ber ] ‚ LE, 


q =0,J<n 
which means from (6) 
(8) in) = su m 
and 
(9) qua) m ea Pru @). 


To determine the remaining constant c,,,, we note that by (7) and (9) 


(lx?) z Pe P) — & (1x) POR) BAZA =e Be 


AS 


Letting x = 1, then from the hypergeometric series definition of the Jacobi polynomials 


we have P* (1) = un. which gives 
Ger (a), i 
— en 
and so 
(10) ah 


Combining (5), (6), (9), and (10), we have 

l x, x a=1, p- — = 
(11) p ED. PM) (lx) (14x) ]=-2 E á x) (1—x)* (+ il 
To obtain the second order ordinary differential equation, change 5 — n—1, « — a4 1, 
and fj — fi4-1 in (11) to give 


u AAA MEA A A 


dx n1 


Then by (3), 
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1 2 7 [o 4 Ax A x 
A | E m je Do) (1—x) H et | = —2? n P. LES (1—x) (+Ô, 


NU P" satisfies 


dr LE u E |== (n+a+ß+1) (1-2) (1239 y. 


Using the product rule to expand the left-hand side, we have 


(1—x) d |-(«+1) (1-+x) y’ a (B+1) Usd JE (1—x?) | 
= —n (nta+fp+1) (1—x)* (14x) y 


which, when simplified, becomes the second order ordinary differential equation for the 


facopi polynomials y= P(x) : 
(12) (1—x^) y^ t [ (f—a) — (a--2) x ]y' - nno fi-1) y 2 0. 


The reader is invited to compare this result with the second order ordinary differential 
equation for the Chebyshev polynomials (T.(x)) given by Equation (5) in Chapter I. 


By iterating (11) k. times, we obtain 


k 
A [ P P6) Ur Y =D AP A A 


dx 


Setting a >a+k, and 8 — +k gives 


x A CF i x a 
(1—x) (1+.x)? Po) == 22 [(1—x ) +k fee. s ps +k, — |. 
2 (14-1), dx 


or equivalently in terms ofthe weight function, 


E a 
a (13-1), dx" 





w(x ;o, f) P P) — [ w(x sath, B+k) Pet Pr] 


the general Rodrigues’ formula for the Jacobi polynomials. Letting 7 =0 and noting 
Py (x) = 1 gives 





X A En) d® X 
E OS 7 [ (15:29 ** (14597 ** ]. 
CUN. AN 
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OI 


AA 


(x, f) 
wre E E 
` V 


| w(x; a+k, B+k) ], 


the classical Rodrigues’ formula for the Jacobi polynomials. 
Finally, we can use these ideas to obtain the value of [A^]. By direct com- 


putation, 


b ime = [ (1—x)” (1+x) dx 
EI 


2 oe pm 
[ (z+ B +2) 


We will use this below. Combining (8), (10), and (4) when m =n yields 


a, fi A(n+1) "L 
D a ie Dom f 2] de 


Making the changes n > n—l, œa > a+l,and £ — fi--1, this may be rewritten as 


Ele - Her | j: 


Iterating this relation n. times produces 


[ít e H (n+a+ß+1), P 


4” n! 
_ T(2n+e+D-+1) oon vay rco TIS eee 
2^" n! T(n+a+8+1) ido pana) 


A a 
(2n+a+8+1) an! Dinta+ß+1) 


3. Generating Function 
As was mentioned in Chapter II], finding a generating function for a class of 
orthogonal polynomials can be a challenging task. Fortunately, a generating function 
for the Jacobi polynomials can be found by mimicking a technique attributed to [fermite 
in his work with the Legendre polynomials. (See Chapter TL Section too) 


We seek a generating function 


o per (x) r” 


As 


(13) DM 


I 
c 


n 


femeonstants c,, 7—0, I, 2,..., with 


| pe PC x*(1— x (L4 xo dx =0, k <n. 


To begin, we multiply (13) by x‘ (1 — x): (1 + x)? and integrate from -1 to 1 to obtain 
k 
E! k 1 ß (a, B) X ß n 
a4) x° F(r,x) (1 — x) (1 + x) dx=) c, [ jor (x) x“ (I —x) (I + x) dx |” 
E n=0 E 


Note that the summation of the right-hand side 1s from 0 to & because the orthogonality 
property makes each term zero for n >k. We observe that the unknowns to be deter- 
mined in (14) are F (r,x) and c,. When one is given, the other can be found, so we first 
consider the left-hand side of (14) 


| x F (r,x) (I —x) (1 + x)? ae 
= 


l Ca 
Setting /1-2xr+r? = 1—ry and substituting for x 2 y + = r vields 


l (1-y°) r |k (1-7) r U | (1—y°) r | | (1-y°) r | 
Ze ish, ee) or E l — = |=) dy. 
T E 5 | Va y > Be ee: (1779) 4 


Factoring out the terms (1—})* and (I+y)?, we have 


1 ]—y? k (127 
(15) 7 
: re A 


We would now like to make a judicious choice of / in order to facilitate calcu- 








lations. Following the lead of the Legendre polynomials (Chapter III, Section E.3) 


suppose F is such that 
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Is mTC a 
(16) (ne Vf un | LE d ge 0-9) = j: 


then the integral in (15) becomes 


1 l— : Y 
J E cM Y ey dy. 
1 


We note that the integrand is a polynomial of degree k in the variable r as 1s 


the right-hand side of (14). From (16), we get 


I=y* ar use ELS 
rr y EA 7 )-a-sr p 237 | == 2 i 


which becomes, via l—2xr+r = (l—ry}, the generating function for the Jacobi 


polynomials: 
F (rix) 2 2**8 (1=2xr+r 1? [1-7 (1-2xr-r?) "T? [1+r+(1—2xr+r "T. 


Using this generating function in (13) gives 


2 Cr p "= 2 (1—2xr+r [17 re (1-2xr- 7)! ^T? [Lr (1-2xr- 72)! ^ T? 


(a4- 1), 
no c 





n=0 


from which we can determine c, by setting x = 1 and recalling that P® (1) = 
Thus, 


Gr Da. py Ip) oq 


n! 





As 


Cn 
0 


n 


which bv the Binomial Theorem becomes 








and so c, 2 1 for all n. (Ref. 8] 
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B. SPECIAL AND LIMITING CASES 

With the structure of the Jacobi polynomials established, we turn now to the role 
of the parameters c and f. 

]. Special Cases 

For certain choices of the parameters a and fi, we find that the classes previously 

examined and several new classes are produced as special cases of the Jacobi class. 
These subclasses inherit the structure of the parent class which often provides a direct 
way to establish specific properties (1.e., polynomial nature, orthogonality, etc.) 


Table 1 provides the choice of parameters for selected classes. 


Table 1. PARAMETERS FOR SELECTED CLASSES 


Jx obi a>-l, fi25-1 


Chebyshev, E 
First Kind 

Chebyshev, 

Second Kind 











Legendre 






2. Limiting Cases 


In this section, we briefly examine the Laguerre polynomials and the Hermite 
polynomials. Using the hypergeometric series definition of the Jacobi polynomials, we 
show how the Laguerre class is a limiting case of the Jacobi class. 

a. Laguerre Polynomials 

The Laguerre polynomials, L(x), are defined in terms of hypergeometric 


series 


na —H 
i i-us] 
n = pe) | i 


The first relationship to establish is 


ol 


a . x; D 
LS (x) = Ze ps 2) (1- EN ) 


To show this, we begin with the hypergeometric series representation for the Jacobi 


polynomials 


a naa —n,n+ta+ß+l l-y 
Pe Py) =( ar a 


n a+ | 


When y = I— Es , we have 


b 
n+a —n, ntatp+l 
pe P (1-2- )-( m p = 
b n a+] A 
Writing out the power series, we obtain 


ph (1-2) - (n9) 
Í p a k 


(—n), (n+a+0+1), x 


s 


= (a+1), k! pe 
7 ( n+a ) y (=ni (nF 
n /£4 (a+1) k: he 
In the limit as $ > oo, the ratio 
(n+a+ß+1), 
D =l 
and so 
, 2x n+0 ` (—), L. ) 
la rem (1-28 )=( J A A 
Boo i ß H = (+1), k! C ) 


completing the proof. [Ref. 13: p. 103] 
To obtain the orthogonality relation for the Laguerre polynomials, we start 


with the orthogonality relation for the Jacobi polynomials 


l 
J PEA DP AY Ue qy—0, msn. 
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2x 


p 


p 2 2 2x Ya B2 
pe k = Joe m1 x F )| —dx=0, men. 
| A 5 5 f f x MEN 


Passing the constants outside of the integral and dividing through by them leaves us with 


B 2 p X 
pe "hi- 3E) rt "hi- 3E s (1- Y axo. m = hn. 
| m pj i ß 


Taking the limit as 0 — oo, we obtain the orthogonality relation for the Laguerre 


Letting y — 1— gives 


polynomials 


(18) | LE) LE) x'e"^dx20, msn. 
0 


b. Hermite Polynomials 


The Hermite polynomials, H (x), are defined 


m E 


(19) IL(x) — (2x Ah] onu | 


In a fashion similar to that for the Laguerre polynomials, the Hermite polynomials are 


a limiting case of the Jacobi polynomials via the Gegenbauer polynomials. Specifically, 


IT, (x) = n! lim e ME 


. >00 


which allows a derivation of the orthogonality relation 


D^ : 
(20) | EMEN EE s 


from that of the Jacobi class. (Ref. 13: p. 107] 


C. DISCRETE EXTENSIONS 
We turn now to orthogonal polynomial classes which use the discrete inner product 


introduced in Chapter II, Section A.2 
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N 
<f,g>= à Fe. 
X=0 


Here instead of a continuum the support of the weight function is concentrated of a fi- 
nite set of discrete mass points (0, 1,...,N}. The polynomial classes of particular in- 
terest are the Hahn, dual Hahn, and the Racah polynonuals. 
L. Hahn Polynomials 
The Hahn polynomials - actually discovered by Chebyshev - were independently 
realized by physicists working in angular momentum theory via 3-j, or Clebsch-Gordon 


coefficients [Ref. 8]. We define this class by the gencralized hypergeometric series 


—n, —x, nta i41 
(21) Or, A] : | 


a+l,—N 


for a 2 —1, f» —1, where N is a positive integer and 5 —0, 1, ..., N. From the power 


series 


n 


= —x), (0 l 


k= 0 
we note that the variable x does not appear where we have come to expect. Since 


(—x), 2 (7x) (—x1) ... (-x+k—|) 
= (~1)* (aD)... (KEN), 


we conclude that Q(x; a«a, B, N) 1s a polynomial of degree n in the variable x. Because 
O<n<VN, this set of N+1 orthogonal polynomials is finite for fixed a and f. [Ref. 14: 
Das) 


Ihe Hahn polynomials satisfy the discrete orthogonality relation 


N 

(22) Dy Quixia, f. N)O (x a p ON (X Cn ern 
x=0 

where the weight function is given by a Aypergeometric distribution 
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ee E DS E. 





3 E N—x 
(23) u (a+1), (f -- 1) x 
[p (l 


We note that the weight function may also be written 


(B+1)y (a+1), (EN ) 


E wm By Ne UENDE 


We introduce the first forward difference operator acting on x 
A f(x) = fitx<+1) — f,(x) 


as a discrete analogue to differentiation. Since 


2 Af 0 m f o1 — 40), 
[=Ü 


(a discrete analogue of the Fundamental Theorem of Integral Calculus) the first forward 
difference operator 1s a discrete inverse of the summation operator. This difference op- 
erator is used in the Rodrigues’ formula for the Hahn polynomials, which can also be 
written in terms of the weight function w(x;a, f, N) with shifted parameters as was 
done for the Jacobi polynomials. 

Two limiting cases of the Hahn polynomials are of particular interest. In the 
first case, replace x by Nx in the interval of orthogonality. This in effect places the 
support of the weight function on the ( N+1) equally-spaced points { 0, I/NV, 2/N,..., 1} 
in [0,1]. As N — oo, the set of discrete mass points tends towards the full interval 
[0,1] ; we expect that this structure is reflected in the discrete orthogonality relation (22 
becoming an integral orthogonality with respect to a continuous weight function on that 


interval. Rewriting the weight function in (23) as 


D(x-1-«4) TOUN—=x+1+f) 


(x: Se 
jubeo mu) (x1) [( N—x4-1) 


and using a consequence of Stirling's Formula [Ref. 15: p. 257] 
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r 
m s ea 
m S 


we see that 


P(a+1) F(+1) lim N77 w(Nx;a, fl, N) = x% (1—x)%, 


which we recognize as the continuous beta distribution for the Jacobi polynomials nor- 
malized on [0,1 ]. 

This is indeed the case, and can be easily verified by a direct computation on the 
hypergeometric series definitions (21) ol Q,(x;a, ß, N) and (1) of PÜ”(«x) given in Sec- 
HOMA T E: 


—1— pe Pl. 


l MO V 
Jim Onl: xX, a, B, i )- zm 


Thus the Hahn polvnomials may be viewed as a discrete analogue and generalization of 
the Jacobi polynomials. [Ref. 14: p. 36] 
The second limiting case gives rise to an interesting class of polynomials which 
has applications in coding theory [Ref. 16]. 
a. Krawtchouk Polynomials 
For0<p<l,leta=pr, ß=(l-p) tin (21), then take the limit as ¢ tends 


to infinity to obtain the Krawtchouk polynomials, 1.€., 


lim O,(x;pt,(1l-p)t,N)= 


25 puer. CIE 
( ) 2 | EM , p IE nM SD. | 


[Ref. 14: p. 38]. A similar limiting argument applied to the weight function shows that 


the Krawtchouk polynomials are orthogonal with respect to a binomial distribution: 


(26) Ste: PEO o Jr (1 -p) * 290, m x n. 


x=() 


This class of polynomials possesses an inherent symmetry, the structure of 


which can be generalized to form other orthogonal polynomial classes. To this end, we 
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turn now to the characteristic of duality. For suitably defined functions p(x) and v(x), two 


classes of orthogonal polynomials { p,(u(x)) }, and | q,(v(x)) } are said to be dual if 
px(utn)) = gy (v(x)). 


That is, interchanging the roles of the degree n and the discrete variable x in one class 
produces the other. The Krawtchouk polynomials provide an example of an orthogonal 


class that is self-dual, 1.e., 
K(x; p, N)=K(n;p, N). 


This is clear from the hypergeometric series definition (25) given above. It therefore may 
seem reasonable to suspect that there exists a class of orthogonal polynomials dual to 
the Hahn polynomials. Such a class does in fact exist, and it is this dual class which we 
Mexi examine. [Ref. 17: p. 657] 

2. Dual Hahn Polynomials 


The dual Hahn polynomials, R( A(x) ; v, ó, N), are defined 


—n, —x, X--y T4] 
(27) (y o, N | 


v+], — N ' 


where ¿(x)= x(x+y+ó+1) and n=0,1,..,N. These polynomials satisfy the 


orthogonality rclation 


N 
(28) p O EN RA one, N) Omen 
x0 
where 
| (yt5+1), ((y+5+3)/2), (Vt), (CN ). (-1)* 
A Ld em 
(29) (1), ((y+-0+1)/2), (O+1), (y+0+.NV42), 


(„Hör (vH. N.  2xty+d4+l 


= 2200 XI a > 
p ron DON 7,02), y--à4-1 


Note that A,( 4(x) ; v, , N) is a polynomial of degree n. in the “variable” A(x). The rca- 


son can be scen directly from the hypergeometric series 


Du 


n 


Ral Mx); y,0, N)= > 


k=0 


(—n), (—x), (etyt+d+1), 
(7-1 CN Ji K 


By writing the terms (—x), (x+y+0+1), in product form, we obtain 


(—x+/) IM (O EIN 
/=0 j=0 
which becomes 
k—1 
| [erp (r+y+ó+14)) 
j=0 


Multiplying these factors as binomials, we have 


k-} 


| | [ (x) ety tot) $j (Hx) ti (et yt 541) +] 
J=0 


which simplifies to 


k—1 


| [ —<(x+y+ó+1) +j (y+ó+ 1) + À ]. 
j=0 


Taking only that part which depends on x yields 
A(x) = x(x+y+ó+I) 


as given above. (Ref. 18: p. 48] 
As discussed in the previous section, the discrete classes interact more naturally 
with difference operators than with differentiation. To accomodate the quadratic form 


of A(x), we introduce the divided difference operator 


A fj AG) -— fa 46er D) f Ax) 


P Jnl A(x)) = Al) — A(x-10)-4(x) 


The dual interplay between the dual Ilahn and IH lahn polynomials follows easily 


from their hypergeometric definitions: 
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—x, —n, n--a-- f-4-1 
MIS Boro jÜ; N). 
en | Qx; a, B, N) 


R (A(n);e@e, B, N)— 1A 

This duality is also reflected in the recurrence relations and difference equations of the 

two classes. By interchanging the roles of n with x, « with y, and with ô, the recur- 

rence relation for one class leads to or can be extracted from the difference equation for 
the other class. [Ref. 14: p. 37] 

3. Racah Polynomials 

Both the Hahn and dual Hahn polynomials can be unified as special cases of a 

single larger class, the standard notation of which is similar to that of the dual Hahn 

class. The Racah polynomials, R( A(x) ; «, B,y, 6), are defined by the hypergeometric 


series 


—n, n+atB+l, —x, x+y+6+1 | 


30 R,( A(x); 0, B, y, 0) z 4E. 
E nh AX) 5 0 Ps 70) = 4 | atl, B+d+l, y+! 


where A(x) 2 x(x-4-y--ó--1), 1 2 0, 1, ..., N, and one of a+l, f+6+1, or y+1 equals —N. 


(Physicists understand these objects via 6-} symbols.) 


The orthogonality relation for the Racah polynomials is 


N 
3) Ral Ale) 5a, Boy, 5) Ry AQ) sa, By, 5) w(x 5a, B,y,5)=0, men 


X= 
where 
EN (y+6+1), ((y+043)/2), (a+1), (B+0+1), (y+1),, 
b wars ae Pos = SD) j C HD), 
62 ++) (a+ (B+0+1) (VD Deby +541 
~ (1), (ytd-at+l), (y-B+D, (641), ytd +1 


[Ref. 19: p. 24]. Note that w(x; a, f, y, 6) has a "well-poised" structure. This means that 


the pairwise sum of numerator and denominator parameters is constant, Le., 


(y+0+1) + (1) = y+04+2 

y+ô+3 yto+l 
(28809) (2B) gia 
(a+1) + (y+0—a+1) = y+0+2 
(P+0+1) + (—B+1) = y+0+2 
+1) + (ó+ 1) = y+0+2. 


The same could be said for the weight function of the dual IIahn polvnomials. (To say 


that a hypergeometric series 1s well-poised means that for 


F a), e Ary . 
r-lr b b y Y , 
po Dr 


the parameters a and 6 satisfy the following relation 
a,+l = a,+b, — — ame 


Well-poisedness is an important property of certain summable hypergeometric series.) 
From the hypergeometric series definition (30) of the Racah polynomials, we 
find that when the roles of x and n, « and y, and f and ö arcall ınterchanged, the series 
is unchanged. Thus, like the Krawtchouk polynomials, the Racah polynomials are self- 
dual. 
To recover the Hahn polynomials as a limiting case of the Racah polynomials, 
let y--1 2 —N and à — co. Thus by formulas (30) and (21), 


lım 47; 


Ô> 


—n, n+a+ß+1, —x, x+y+ô+1 T E —n, nta+f+l, —x 
atl, B+d+l, V E pee a+l, —N 


and then by formulas (32) and (24), 





l 
dim w(x; a, $, y, ô) = TR w(x;a, B, N). 
Likewise, by letting «+1 = —N and f > oo, we obtain the dual Hahn polynomials. 


Figures 5 and 6 together with Tables 2 and 3 provide a hierarchy of classes dis- 
cussed in this chapter. Tables 4-15 summarize information about selected classes in this 
hierarchy [Refs, [2098197 20, 21505 9159 2I 
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Racah 


Hahn Dual Hahn 


Jacobi Krawtchouk i Meixner 


N 


Laguerre Charlier 


Hermite 


Figure 5. Hierarchy of Limiting Cases 
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Jacobi 


v 


Gegenbauer 


Chebyshev Chebyshev Legendre 
First Kind Second Kind 





Figure 6. Hierarchy of Special Cases 


Table 2. TABLE OF CLASSES 


Symbol [[ypergeometric Series 


Rel Aa) sa er 
Ry A(x) 5745, N) ee | 
E PE e 


Krawtchouk KAP N) 


Meixner Neo.) 


| 112, (—=n+ 1/2, 1 
X 
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Table 3. TABLE OF CLASSES (CONTINUED) 


M 3 = 
Gegenbauer ln r Bird | 


Asa ° 2 


Chebyshev, 
First Kind 


Chebyshev, 
Second Kind 


Legendre 
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Table 4. RACAH POLYNOMIALS 
Symbol: R(AÀ(x);a, f), y, ó) 


us v l bor a positive integer, x — 0, T, ... , N. 


Weight: 
_ (yto+l) ((y+o+3)2)  (z+1), (6+ó+ 1). (y+ 1), 


vert Pr) = ay (eie. NEN (AD 


Norm: 
N 
2, [RGGQ ias Boy) P wGrs es $, y, 0) 
ss 
n! (n--&-- fi--1), ( B+1), («-d4+1), (a+ B—y+1), 
= 1 f E E EA EN TTA A E e 
(LP SIDE Eo yd 
where if 
+l = — N, 
y ++ Cy 
(y—f-- 1) (6-1), ` 
or if 
ß+ö+1=-N, 
" (y+0+2) y [ó-0) y 
en EIE 
or if 
y+] = — N, 
l M (—ô)y (&+8+2)y 


G js | 
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Table 5. RACAH POLYNOMIALS (CONTINUED) 
Hypergeometric Series: 


; =i, +B+1, —x, x+y+6+1 
Ral A(x); a, B, y, D =R i ç. i E 4| 


where A(x) = x(x+y+ö+l), n=0,1,.., N, 
and one ofa+l, f+ö+l, or y+lequals —N. 


Recurrence Relation: 


A(x) Ral A(x) 50, By, 6) =aln)[ Rag A a, B, y 5) ] 
—c(n) [ RC A(x) @, B, y, 9) Ry, (460) 5o, B, y, 6) ] 


where 
(i s Mob Esp s n n si 
n (2a+a+6+1) (2n+a+f+2) 
A IN) 


(2Qnta+fB+1) (2n+a+f) 


Difference Equation: 


An) RA Ala): a, B,y, 5) 2 AGO | RC AG 5o By, JRC A(n) 5 a, B, y, 5) | 
CR) RL A), BO) RI AG a Boo 


where 
pm (x+y+0+1) (x+y+1) (x+Ó0+8+1) (x+a+1) 
7 (2x+y+6+1) (2x+y+6+2) 
x (x+6) (x+y+d—a) (x+y—f) 
C(x) = ————[——,—,ÀOO "““S.HE———Ə——əs 


Qx+y+0+1) (Qx+y+0) 
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Table 6. DUAL HAHN POLYNOMIALS 
Symbol: R,( A(x); y,6, N) 


Interval: For Y a positive integer, x =0, 1,..., N. 


Weight: 
aer, yy e LEPE (OD + DN) E 
CUPIT (eG AA 
Norm: 


N 
4 N— El 
2 LRQ 93, N) wi 3, | (1 Pus) | 


x=0 


Hypergeometric Series: 


—n, —=x.x+y+ó+ 1 
R.(À(x) ; y, 6, N] " n b 4| 


where A(x) = x (x +y+ó+]) 


Recurrence Relation: 
—A(x) Rf A(x) S O, N) = D(n) [ Bl A(x) yt O, REN A(x) > Vy Ó, N) ] 
— D(n) [ K A(x) Ys Ó, N Er A(x) Tr Ó, N) ] 
where 


B(n) = ( N=n) (y+142) 
D(n) = n( N+1+6—n) 


Difference Equation: 
—aR,( A(x); y, 6, N) = d(x) [ Ry B ON S O yo T 
E [ RC A A p o y Je Ao) eo N) | 


where 
b(x) (x+y+ö+l)(x+y+l)( N—x) 
x) = ———— — ue—— 
2xty+ö+l) dxr+ty+ö+2) 
x (x+0) (x+y+04+N+1) 
MINE == 


Qx+y+0) Qx+y+0+1) 
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Table 7 HAHN POLYNOMIALS 
Symbol: Clr a, ey ) 


Interval: For N a positive integer, x 2 O0, 1, ... , iV. 
Weight: 


"y | f(x*taN ( N—-xfN C 
w(x ; a, p. N ) -( X ) ( Nx )- Tp (D, —-N-f). (-N-f), 


Norm: 


? [Qs f, N)y w(x;a, f, N) 
x=0 


lier di E" u) 
oe Gef N n 


~ 2n+a+ß+! (\ 
H 

[ (a+ 1) [ (z+ + 1) L(n4- B-1) FG 1) 

[(f4-1) I(n+a+1) L(2-a4- f-1) 


X 


Hypergeometric Series: For a, f» —l and 
N 


zur N 
| 
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Table 8. HAHN POLYNOMIALS (CONTINUED) 
ISecurrence Relation: 
—x Q(x;a, f, N) 2 b(n) | Orta. pdt) Or Ger om, N) ] 
un) [On sa, B, N) - Quaeso, f. N)] 


where 
b(n) = (n+a+B+1) (n+a+1) ( N—n) 
P Onta+f+D) Qnta+f+2) 
Ay n (n+ß) (n+a+ß+N+]1) 


(Qntat+f) (2atatp+1) 


Difference Equation: 


—/(n) O,(x;0, B, N) = BG)| Q,(xo-3;o, B, N) —Qu(x;o, f, N)) 
—D(x)  Q,(x;a, B, N)  Q,(x—l;o, fl, N)] 


where 


B(x) 2 ( N—x) (a4-14-x) 
D(x)=x( N+1+ß-x) 
An)=n(n+a+ß+l) 


Rodrigues’ Formula: 


De A 


u ("*P) y | pud [e | 
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Table 9. JACOBI POLYNOMIALS 
Symbol: P(x) 





Interval: [ —1,1] 


Weight: (1 — x)* (1 + x) 


(a+ 1), 
n! 


Standardization: P@ (1) = (" M z) = 


Norm: 


ee 


1 
(a, B) 2 Afs ß Z 
| PAULI AS enn Porte Bad) 
A 


E E 2 E 
Hypergeometric Series: P% P(x) _ (at, Al n, nta+ß+l, 1-x | 


n! a+ l i 2 
Recurrence Relation: 2 (n+1) (n+a+ß+1) (2nta+ß) P9:P (x) 


= (2nta+B+l) [ (2n+a+f) (Qnta+f+2) xta7—p* | PO P(x) 
—2 (n--a) (n--B) (2n--a-- f--2) PU P) 


Differential Equation: 


(=x?) y” +[ B—a— (at f+2) x] y’ +n (ntatf+ly=0, y= pe PY x) 


Rodrigues’ Formula: 


n 


dx” 





Y pe PG) ae 1)” (E ee (14x)? a m De 2: | 


V (1-2x1+1°) 


Explicit Expression: P^ P) 27 Y ELS (s (x—1)7* (x+1)* 


k= 


Generating Function: 

= a+ p 

5 pe Pix na -— (1-7 | — 2x1 + t? y CEN ] — 2.xt + t? 73 
nz 
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Table 10 LAGUERRE POLYNOMIALS 
Symbol: L(x) 


Interval: [0, co) 
Weight: ze”, a>—I 


il: 


n! 





Standardization: 11) = +. 


Norm: 


xor = I" +1 
| x e x dx = inter) 
0 


H! 
Hypergeometric Series: 
(a) ENTA - mm 
Ln (x) = ( n ) F p x 


Recurrence Relation: 


(n--1) L® (x) =[ (2nt+a+l )—x | L(x) — (n--a) Ta e) 


Differential Equation: x y^ -(a--I—x) y'-ny-0, y L? (x) 


Rodrigues’ Formula: 





a l {" A 
p (x) mL =X ; n Ee € >] 
nX e dx 
Generating Function: 
= ee E) exp( = ) 


n=0 


Explicit Expression: 


(x) = k + | k 
LEO A) T a 


k=0 





"el 


Table 11. HERMITE POLYNOMIALS 
Symbol: TH (x) 


Interval: (—co, co) 
Weight: e-* 
Standardization: F/(1) = 2^ x ... 


Norm: 


| EROS li eU dx Jr 2^ p! 


— OO 


Hypergeometric Series 


H,(x) = (2x)" 5 Fl rg 


Recurrence Relatıon: 
HA AN 


Differential Equation: 
y —2xy *2ny-20, y IL(x) 


Rodrigues’ Formula: 


no MEM x? 
IL(x) = (-1)" e”” pro i 


Generating Function: 


y JA u = 2xt=1? 
n! 5 


n=0 


Explicit Expression: 
[n/2] 
s ne2k 
IL(x) e nt». (71) E PLC MNT 


MEE tn 
= kK! (n—2k)! 


d 


Table 12. GEGENBAUER (ULTRASPHERICAL) POLYNOMIALS 
Sabol C(x) (tor P (x), 42-—1/2 


Interval: ( —1,1] 


Weight: (1 — x) 


Standardization: 
x) 
lim ——— = T(x), n=0, 1, 2,... 
SIG OR 
Norm: 


beds o av u 2 424112) Y 
J. LO, (x). ( : ) sa (n+A)n! T(n+2%)) 


Hypergeometric Series: 


DLE, Ra =n,n+24. l-x 
em co | e a 





Recurrence Relation: 


UV Gr a 1) CO) 


Differential Equation: 
(1—x7) y (24-1) x yn (n-22) y 20, y 5 CO) 
Bodteues Formula: CMx)=1, Cf) =21x, 


A CA C ^02)» EG Fin cr 


dx 


Generating Function: 


2? CPG) 1" — (122x187 


n=0 
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Table 13. CHEDYSHEV POLYNOMIALS OF THE FIRST KIND 
Symbol: 7,(x) 


Interval: [ —1,1] 
Weight: (1 — x?) 12 
Standardization: DS] 


Norm: 
| 2 0 
2 23-112, _ jruj2,nF 
| [5,6)] (1 — x^) "dx - S E 
=) 


Hypergeometric Series: 





—n, ] —: 
te | 


Recurrence Relation: T,,,(x) = 2xT,(x) — T,_,(x) 
Differential Equation: (l-—x) y°"—xy'+rn74y=0, y=T,(x) 


Rodrigues’ Formula: 


Gy oa 


= 2er 2\n—1/2 
T,(x) a r (n L 1/2) du^ (1 E ) | 


Generating Function: 
Int eat 
l— = 


Explicit Expression: 
T,(x) = cos(n0) with x2cos0, 00m 


[n/2] 
T,(x) = £ E B 3 C 


k=0 


74 


Table 14. CHEBYSHEV POLYNOMIALS OF THE SECOND KIND 
Symbol: U (x) 


Interval: [ —1,1] 
Weight: (1 — x2)1⁄2 
Standardization: U(l)=n+1 


Norm: 


1 
| [U,(x)]? (1 — — = 


Hypergeometric Series: 


U 9) e ena 7t | 


Recurrence Relation: U,,(x) 2 2xU (x) — U, (x) 
Differential Equation: (1 — x2) y^ — 3xy' +n (42) y =90, y= U(x) 
Rodrigues’ Formula: 


(-1)"(n +1) /x q" 


L cae 
(Lax?) "7 24! P(n4-3/2) de 


n(x) = 





a ae 


Generating Function: 


> U l = — SA ESI 
=0 1 -2xt--t 


Explicit Expression: 
sin(n 4- 1)0 


us sin 0 


with x=cos0, 0 < D «or 
[1/2] 
! CE n-2k 
UG) = > (-1) "Un — 28) (2x) 


k=0 u 


> 





Table 15. LEGENDRE (SPHERICAL) POLYNOMIALS 
Symbol Pl) 


Interval: [ —1,1] 
Weight: 1 
Standardization: P(1)= 1 


2 


Norm: í [P(x] dx = TES 


Hypereeometric Series: (x) =F EN | 0 


Recurrence Relation: (241) P,.,(x) = (Qn+1) x PL) — n P,_¡(x) 
Differential Equation: (1 — x7) y" — 2xy' c n(n-1)py 30, y= P(x) 


_1Y 
Rodrigues’ Formula: P(x) = - D < [ (1—x2)'] 


Generating Function: 


5 P (x) i” = (1-24)? , Jere DP i 


n=0 


[n/2] E 
Explicit Expression sr o) E Y (-D n E ee 


n 
2 k=0 n 
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V. APPLICATIONS 


Orthogonal polynomials and special functions in general have been studied exten- 
sively in mathematics and other fields since the eighteenth century. Presented below are 
a few of the traditional applications of selected classes of orthogonal polynomials. Our 


first few applications come from numerical analysis. 


A. ECONOMIZATION OF POWER SERIES 

Economization of power series is a technique used to reduce the degree of a 
polynomial approximation to a given function. 

The maximum norm (or L^ — norm) for a continuous function on a compact inter- 


val [a,b] is defined as 


= ma: XJ. 
I Heo = max Io)! 


This norm is not induced by an inner product, but nevertheless has many uses in applied 
mathematics, including numerical analysis. 

The minimax property of the Chebyshev polynomials states that of all »^ degree 
monic polynomials (1.e., leading coefficient 1), 2'^^7;(x) has the smallest maximum norm 
on [ —1,1] [Ref. 23: p. 106]. The justification for this statement is deferred until Section 
C. Fence the best approximation in the maximum norm to the function x^ on [ —1,1] 
bv a function of lower degree is f(x)2."—2'T(x) So, given a function and a 
polynomial approximation to that function (e.g., from a Taylor series expansion), suc- 
cessively replace the highest powers x^ with f/(x) to obtain a polynomial approximation 
of lower degree (Ref. 23: p. 125]. 


Example 1: Let f(x) 2 sin x. The Maclaurin series for this function 1s 


oO 


k 

| (m) 2k+1 
5 l 
2, (2k + 1)! 


and this series is convergent for xe R. IF truncated after the x term, the polynomial 


approximation for the function is 


3 5 

; x x. 
SIn x x= E en 
3! 3! 
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with a maximum error of 0.0002 for xe[ —1,1]. We use the fifth degree Chebyshev 


polynomial 7;(.x) = 16x5 — 20.9 + 5x to obtain 














x 2 5x” zada 
SIS 1920 
with an error not exceeding 
max | 75(x) | ] 
— 1950 = 000 * 0.00052 
in [ —1,1]. Thus the approximation 
3 3 
SEN CI e O NIE EDI ee 75 
~ 3! 480 1920 1920 480 


has an error whose magnitude in [ —1,1] does not exceed 
0.00052 + 0.0002 = 0.00072. 


Compare this with the maximum error of 0.00833 for the Maclaurin series which is 
truncated after the x? term. For a cubic polynomial approximation of sin x, the “econ- 
omized” polynomial has a maximum error that is significantly smaller (less than one 
tenth) than that of the truncated Maclaurin series. 

The next three applications illustrate the usefulness of the zeros of orthogonal 
polynomials. The first two come from numerical analysis, the third from a problem in 
electrostatics. We begin with a preliminary discussion of a fundamental technique from 


numerical analysis. 


B. POLYNOMIAL INTERPOLATION 

Polynomial interpolation is a method of approximating a given function with a 
polynomial that matches (interpolates ) the function at specified points (called nodes or 
abscissae) x,,...,X, [Ref. 24: p. 497]. Given a function f(x) and 5 distinct nodes in a 
compact interval [a,b], there is a unique polynomial of degree (1 — 1) that passes 
through the points f 0 le eee 


l'or each /— 1,2, ...,7t, define a polynomial of degree (n — 1) by 
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n 
| | X— Xg 
mi X) = ( X), ) 
k=] 
ki 


Clearly, for each k =1,2,...,n, 2,,(x,) =0 if A#i and xz, (x,) = l, that Is, x=, (x,) = 0,4. 

(We may thus equivalently express 7,,,(x) = J,(x)/((x—~,) 1’(x,)), where L(x) 2 [] (x—x) is 
k=1 

the unique monic polynomial of degree m that vanishes simply at each node x, 

Heel 2, os. ) 


The Lagrange interpolating polynomial is now given by 


(x) = > f(x) mpl): 


i=!) 


Clearly, L/(x,) 2 f (x,) for k 2 1,2,..., n, and uniqueness of the degree (n — 1) interpolat- 


ing polynomial is guaranteed by the Fundamental Theorem of Algebra. 


C. OPTIMAL NODES | 

In this section we address the issue of estimating the maximum size of the interpo- 
lating error || f— L/ ||. Assuming fis suitably differentiable in [a,b], it can be shown 
[Ref. 25: p. 188] that there exists a value K, (which depends on fin [a,b]) such that for 
anv x in [a,b], 


K 
| ix) | 





Ia Enon = 


muse) - [T (x — x). (Note again that Z(x) 1s a monic polynomual of degree n which 
vanishes at the nodes.) This implies that || / — Lf || is minimized by making the optimal 
choice ol nodes x,,..., x, in [a,b] which minimizes | 7, |... Surprisingly perhaps, this hap- 
pens precisely at the zeros of the Chebyshev polynomials T;(x), scaled to the interval 
[a,b]. We sketch the reasons below. 

For simplicity, we take our interval of interest to be [ —1,1] instead of [a,b] without 


loss of generality. The transformation 


u D 
“A - )«( a j 








or cquivalently, 
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x—a 
ds 2( b—a i E 
Is a one-to-one continuous mapping between the intervals [ —1,1] and [4,5]. 
Now let x, be such that 7,(x) =0, ¿=1,2,..., ¿(see Chapter l, Section Aro venues 
ES X T(x). From the recurrence relation (Equation (4) in Chapter 1), T,(x) has a 
leading coefficient of 277! ; hence this /(x) is monic. Moreover, from the definition 


T.( cos 0) 2 cos n0, it follows that 











; l 
(1) lO) T CD 5 
COS a : : : 
where y, 2 y,, =—— ,i=0,1,...,n. Now suppose m,(x) is another monic polynomial 
ol degree 7 such char 
l 
(2) | Pt leo < I fa Ieo m er: 


Combining (1) and (2) we see that we must have for i 2 0, I, ... n, 


m,y) < hly) = = if i is even, and 


(3) E 
my) > hy) = T if ¢ 1s odd. 


Thus, the polynomial px) 2 m,(x) — L (x) has degree at most (» — 1) (since it 1s a differ- 
ence of two monic polynomials of degree nj. with at least n. zeros, one in each interval 
(Y. Ya) Qian- l1, by (3). This contradicts the Fundamental Theorem of Algebra, 


and so no such polynomial m,(x) satisfying (2) exists. Ilence the choice of nodes x, de- 
l 


> n-i 


polynomials /(x), and therefore over all possible choices of interpolating nodes x,. 


termined by /(x) T(x) minimizes S max | I(x) | over all possible monic 

It should be emphasized that this Chebyshev interpolation allows an a priori error 
bound for all x, but is not always best possible for every x using other interpolation 
schemes. (For example, if equally spaced nodes xi- , x, are used, then trivrally there 
ror /(x) — L/(x) 2 0, even if x 2 x, is not a Chebyshev zero.) Clicbyshev zeros are onur 
mal when one has freedom in the choice of nodes. Tor details regarding the practicality 


of Chebyshev interpolation and some asymptotic results, see [Ref. 24]. 
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D. GAUSSIAN QUADRATURE 

Quadrature formulas are used in that area of numerical analysis concerned with the 
approximate integration of a function f(x) against a weight function w(x) 7 0 on an in- 
terval (a,b), when the explicit evaluation is intractable. 

An interpolatory quadrature 1s such a rule that uses interpolating polynomials, such 


as Lagrange polvnomials: 


y b 
| f (x) w(x) dx «| L(x) w(x) dx, 


where L/(x) 2 f(x) n(x), as described above. This can be rewritten as 
iz] 


b n 
(1) | reo wao S) w 
a i=! 


where the weights w, are given by 


(2) = | Ale) wi) dr 

a 
Thus, interpolatory quadrature 1s basically a weighted sum of the function values f(x,) 
at the nodes x,,i=1,2,...,n, as are numerical integration recipes such as Simpson's 
Rule and the Trapezoidal Rule. 

Bornspecified nodes x, ... ,x, e 14,5], then weights w,, ... , w, computed in (2) for the 
quadrature (1) will be exact for polynomials fe P,_,[a,b], but we can do better. [Ref. 
2555-236] 

In Gaussian quadrature, we ask for the location of the » nodes x, ... , x, as well as 
the n weights w,, ... , w, in order for the quadrature rule (1) to be exact for polynomials 
fe P,, ¡[a,b]. At first glance, this seems to be an extremely complicated computational 
problem, but the solution falls out simply when the theory of orthogonal polynomials 
is applied. 

Let ! p,(x) 1%, be the class of polynomials orthogonal with respect to the weight 
function w(x) on [a,b], say, by Gram-Schmidt, and let L/(x) be the Lagrange polynomial 
that interpolates fe P, ([a,5] at the zeros of p(x), so that f (x) — L'(x) e P,, [a,b] vyan- 
ISmestar Wile zeros also. Since the zeros of p(x) are real, simple, and lic in (a,b), we have 


the property that 


8l 


yao) zz L(x) =; Pr) q. 


wherc q,., € P,_.[a,b]. Hence, by orthogonality, 


b b b 
| f (x) w(x) dx — | L(x) wld = | p EET AN 


We have thus shown that the resulting n-point Gaussian quadrature rule is exact for 
fe Pala, b]. Moreover, by a theorem of Stieltjes, if f(x) is continuous on a finite inter- 
val [a,b], then 


b 
lim U(x) w(x) dx = [ f (x) wx) dx. 


a 


By applying the Christoffel-Darboux formula (Equation (9) of Chapter III with 
notation from Chapter iII, Section B) and using (1), it 1s possible to derive an alternate 


expression for the weights (2): 


= Ayzı 
hy, Pn (Xj) Boule 


which are referred to as the Christoffel numbers. [Ref. 8 ] 
Thus Gauss-Jacobi, Gauss-Chebyshev, and Gauss-Legendre are the names given to 


Gaussian quadratures involving the weight functions and orthogonal polynomials from 


the Jacobi, Chebyshev, and Legendre classes, respectively. 


E. ELECTROSTATICS 

The zeros of the Jacobi polynomials play an interesting role in a problem of Stieltjes 
concerning electrostatic equilibrium. In this problem, fix “masses” of positive charge « 
and $ at the points x = | and « = —1, respectively. Then place n point masses of positive 
unit charge in the interval (—1,1) so that they are free to move. These interior masses 
arc now subject to a “logarithmic potential”, that is, a repelling force that 1s proportional 
to the logarithm of the distance separating them. The problem is to determine the dis- 
tribution of the point masses x, 7 = 1, 222m wlenwlie systems Ee GRUT LES ETCURNE 


Mathematically, this is equivalent to maximizing the force function 


n n 
G) Fly, %,) =a) log(I—x) +B > log(ttx)+ 2 log | xs, 
vs] i=1 


1<i<j<n 


The logarithmic terms give the restrictions x, € —1, x, € 1l, and x, x x, for iJ. This 


otherwise continuous function gives the equilibrium points by setting 


(4) == E et 





and solving for x,,...,x,. To solve this system of n nonlinear equations in n unknowns, 


Stieltjes introduced the polynomial 


and reduced (4) to 








which becomes 
(1—x^)p," (x) -- [ 28 — 2« — Qa -- 28) x, ] p, (x) = 0 


Eu [. 2....,H. 


Since the polynomial 
(1- x?) p(x) - E25 — 2a — (2a + 28) x ] p (2) 


is of degree at most n and vanishes at x = x, i= 1,2,...,n, it can be set equal to a scalar 


multiple /p,(x) which also vanishes at these points. [lence 
(1— x^) p,"G) -- E28 — 2« — Qa + 28) x ] p,'(x) — jp, (x) = 0. 


Attempting a power series solution to this second-order difTerential equation leads to the 


observation that polynoniual (i.e., terminating) solutions exist if and only if 


Ze oS) 
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Rearranging terms yields 


(1 — x°) p. (x) + [ (28 1) — (2a = 1) - (Qa — 1) + EE ) + 2) x | p,'(x) 
tn(n + (2a — 1) + (28 — 1) + 1) p,(x) = 


which is the differential equation for P? "Po. Thus the equilibrium positions of the 
unit charges occur at the zeros of this Jacobi polynomial. [Refs. 8, 13: p. 140] 
The zeros of the Laguerre and Hermite polynomials can be developed as solutions 


to similar electrostatic equilibrium problems. 


r. SPHERICAL HARMONICS 
We investigate another application of orthogonal polynomials to problems in 
mathematical physics. In Cartesian coordinates (x,y,z), the Laplacian operator of a 


function u(x, y, z) 1s defined as 








When converted to spherical coordinates (7, 8, @), this operator acting on a function 


u(r, 0, bp) becomes 





TR 2 l : l 
(5) Vu= E 6 ur), + ENS (sin 0 uə)ə + RS wo | 


where subscripting with r, 0, or $ denotes partial differentiation with respect to that 
variable. A function u is said to be harmonic in a region D if in that region it satisfies 


Laplace's equation: 
Vu — 0. 


In particular, if the boundary of D 1s the unit sphere centered at the origin, then using 
the method of separation of variables, the Legendre polynomials arise naturally as part 
of the solution. 


In separation of variables, we assume that the solution will be of the form 
(6) u(r, 0, $) — f (r) g(0) h(9). 


Substituting (6) into (5) and dividing through by u, we obtain 
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l 


— — h.4— 0. 
h sin°0 = 


ee I . 
(7) T6 f), + EHE (sin 0 gj), + 
The first term, depending only on r, must reduce to a constant which we write as 
v(v-- 1). Substituting v(v+1) into (7) and multiplying through by sin?0, we have 


Lu 0 


v(v--1) sing + LL ( sin 0 g,), + l - hs 0. 


Now we see that the third term, depending only on œ, must also reduce to a constant; 


call it —72. Substituung —7? for this term and simplifying, we obtain 
(8) sin”0 go a + sin 0 cos 0 g, + [ v(v+1) sin20 — m^ ] g ^ 0. 


Dy the change of variable x 2 cos 0, (8) becomes via the Chain Rule 








For problems that are radially symmetric, the @ dependence can be removed by setting 


m = 0, leaving 


d*g 


dg 
dd — 2x —— 4 v(v4-1) g =0. 
x 


dx 





(9) (1-x°) 


When v is a positive integer n, (9) 1s recognized as the differential equation for the 
Legendre polynomials. [Ref. 26: pp. 210-213] 

An alternate approach to this problem (by the method of images) uses Green's 
functions. As motivation, consider the Dirichlet problem for the unit circle in the plane, 
which involves finding a harmonic function u(r, 0) in the unit disk that takes on prc- 


scribed function values f(@) on the boundary r= 1. The solution is given by 


| 27 
ulr, 0) === | f(t) P,(0—0^) q0', 


where 


= 


| —2r cos y+r? 
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is the so-called Poisson kernel for this problem. 


Similarly, the solution to the Dirichlet problem for ihe unit sphere in R? 


Vu=0 
¿= (oni 


can be expressed as 


27 Te 
lr, 0,6) = | |” £0.69 Px) sin 0" at diy 
0 0 


where 

(10) P(y) = e 
[ 1—2r cos y4-r ] 

and 


cos y = cos 0 cos 0' -- sin 0 sin 0' cos($ — d’). 
Writing (10) in terms of simpler functions 


l l 


i + +2r nn | ——————— | 
JJ 1—2r cos yr? zi JJ 1—2r cos yr? 


we note the appearance of the generating function for the Legendre polynomials, with 
x — cos y. Therefore the Legendre polynomials are again part of the solution. [Ref. 27: 
pp. 87-89] 

Laplace’s equation can also be solved in a higher dimensional setting. Let 
X2(x,..,x)e RR". Note then that we may write  x-r/$, where 
r=|xil=Jx?+ .. + x3 and $- P = (č -će R” is a unit vector. A polynomial 
h,{x) 1s said to be homogeneous of degree n if h,(Ax) = 47 A(x). 


The Laplacian operator acting on a function u{x) in p dimensions 1s defined as 


p a 
2 Ju 
M = 


2 . 
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As before, a harmonic function is one which satisfies Laplace's equation Vu — 0. We 
now seek homogeneous harmonic polynomials of degree n in xe R’. It can be shown 


that there are exactly 


2n+p—2 [n+p-3 
N= N, A a ( ) 


n— 1 


linearly independent such solutions, and they can be characterized by Gegenbauer (or 


ultraspherical) polynomials C,(t,), k = 1, 2,..., N. The general solution is given by 
h,(x) = h, (v8) =r" S,(8), 
where the spherical harmonic 


p-2 


—-— ———Á 


N 
SE) = > Ak np cf 2 X. n), 
k=] 


with n, suitably chosen unit vectors. Note that if p=3, then these reduce to the 
Legendre polynomials found earlier. [Ref. 6: pp. 168-183] 


G. GENETICS MODELING 

Karlin and McGregor gave an interesting application of the dual Hahn polynomials 
to a model in genetics. In this continuous time Markoff chain model, the dual Hahn 
polynomials R,( A(x); y, 6, N) arise in the transition probability function for the process. 

The setting for the model assumes N gametes of type a or A and gives a random 
fertilization scheme. The population of either type of gamete is affected by both the 
fertilization process and a mutation process whereby a gamete resulting from a mating 
can mutate into the other type. 

By considering the conditional probabilities for an increase in population size of 
both gamete types, a stochastic process is defined which is a classical birth and death 
process. The transition probability. function for this last process is then cast in terms 
of the dual Hahn polynomials. The interested reader is referred to [Ref. 14] for the de- 


tails. 
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VI. BASIC EXTENSIONS 


A. BASIC HYPERGEOMETRIC SERIES 
[n this chapter, we extend the structure and some of the results in Chapter IV to a 
more general level. This extension is accomplished by introducing a new parameter 


called the base to the hypergeometric series. The base q was used by Heine in a series 


(1-45 (1-4) | (1-45) (1-47) (1-47) (1-477). , 


l l 
E ES (1—9) (14%) (1-95) (1-9 *) 


where c x: 0, —1, —2, ... [Refs. 28,29]. This series converges absolutely for | x | < 1 when 


la| <l bythe Ratio test Since 
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(2) an l-q 


= d, 


we see that the series in (1) tends termwise to the ordinary hypergeometric series as 
q—]. Ihus Heine’s series is called the basic hypergeometric series or the q- 
hypergeometric series. [Ref. 12: p. 3] 

The g-shifted factorial is the basic extension of the shifted factorial introduced in 
Chapter II and is defined 


3 "IER uu 
(3) nn (1—a) (1—aq) .. (1—-aq" 1), "ues 


The ordinary shifted factorial is recovered by applying (2) and (3) in the limit 


GERN: 
mr. — (20% 
elle) 
We also define 
(4) (a;q),, = IN (1—aq*), 
k=0 
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a form we will sec in later results. Since tlie infinite product diverges when both a x0 
and | q | > 1, we will assume | q | < 1 whenever (a; q),, appears unless otherwise stated. 
Hach. 12: p. 3] 

Generalizing the basic hypergeometric series above, we define the ,$, basic 


hypergeometric series (or „p, series) 


Ay V (a, q), iis (a, ; Dn n C) ltr n 

3 Ó; , , x ARE TIO ( (—1) 2 ) x 

e) n VV |- 2 (73 Dn (13 Dn (O55 Dr Á 
1 

where (5) =n(n—1)/2 and g#0 when r>s+1. We require that the parameters 


LM» be such that the denominator factors in each term of the series are nonzero. 


Since 
E O omms m3 2... 


a ‚b, serics terminates if one or more of the numerator parameters is of the form q^" for 


m=0,1,2,... and q xx 0o. When r=s + 1, the expression in (5) simplifies to 


A m q >) (a, 20): Ba aa ar x 
pM Eco. 


Note that in a basic hypergeometric series »' c, x", the ratio c,,,/c, 1s a rational function 
user 12: p. 4] 

Using the g — shifted factorial (3) and (4), we can define basic extensions for many 
of the functions and formulas introduced in earlier chapters. We note that often there 
is more than one way to extend a result; examples will be given shortlv. 


The g-gamma function 1s defined by 


(6) nio. e. Ml 
(1 3 Dow 
Gosper showed that 
(7) lim. (x) 2 F(x) 
q> 1 


<= 1 (Ket. 30; p, 109), Whe stricture of the gamma function extends as well, 


For instance, the formula 
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] — X 
rit =( E re 


can be reduced using (6) to 
I(x+1) = x [ (x). 


With (6) we can define the q-beta function 


which by (7) tends to B(x, y) as q > I>. 
The q-binomial coefficient is defined for integers n and k by 


H G CIT) 
k]a (93 De lGs Mana 


where K 20,1, ...,. For nonintegral « and fj, we define 


B U NEUE 
fi ja m e 
= [ (o1) 


— F4) P1) ` 
The q-binomial theorem is then 


n 


er 
(ab ; q), = 2 H b (aig); qu 
k=0 4 
whene, n0, 12e CIS ESI] 
The next two expressions are basic extensions of the Chu-Vandermonde formula 


(Chapter III, Section F.1) and are both known as the g-Chu- Vandermonde formula: 


5 ss, cq. (c|b ; q), 
( adr == 
2A C i b (c , 4)» 
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p (c 
|". ana EA 


[Ref. 12: p. 11]. These forms can be shown to be equivalent by reversing the order of 
summation. 


Jackson introduced the general form of the q-integral 


| fu) du = | f(0d, — I od. 


where 


[0 ar = atı-a Ita) 
n=0 


The q-integral defines the measure ds which is a natural object for 
q — defined functions (Ref. 12: p. 19]. All of the functions and formulas developed above 


play an important role in generalizing the ordinary orthogonal polynomial classes. 


B. BASIC EXTENSIONS OF ORTHOGONAL POLYNOMIALS 

In this section, we present the q — analogue(s) of selected classes from earlier chap- 
ters. Dy using formulas such as those presented in the previous section together with 
methods based on those outlined in the preceding chapters, it Is possible to derive the 
recurrence relations, difference equations, and Rodrigues formulas as well as many 
other identities satisfied by these q —versions. As mentioned in the previous section, the 
q — extension of a function is not necessarily unique; however the last two classes listed 
are especially important. 

lL. Continuous q —Hermite Polynomials 


a. Definition 


)-» zi. — (q; Dr ín 2098 
q 5 qd) : qux 


k=) 


where x = cos 0. 


21 


b. Orthogonality Relation 


= 2 . 
| H,( COS 0|q) Hp cos 0| q) (e, DL do K lan E 
0 (73 Dac | 
where x= cos d. [Ref. 12: p. 188] 
2. Discrete 4 —Hermite Polynomials 
a. Definition 
[n/2] (OD 


Vu CAECUM 


b. Orthogonality Relation 


n 


[ Hx ;q) H,(x ; q) dy (x) z a) (gig) o e 
=I 


where (x) is a step function with jumps 


| cg Ga 
(Gages 


to 


at the pomtsx = 2, 1=0,12 2. po n on o) 
3. q --Laguerre Polynomials 
a. Definition 


Fora —1, 


C J~“ 
(x) t u Ei? s n+a+1 
L, (x; q) = (q ; q), 191 p »Y,—Xq | 


b. Continuous Orthogonality 


> * dx ea WG soil 
| Di q) DE , q) m a | aC. 1) 0 n 
0 lo , E (50) (usq. | 


c. Discrete Orthogonality 


00 k(x+1) (+ . q) 
Le) k . no k TIE A > no I 
2. nier ie) rer ie) (-c1—44* s" ^ (nod m 


where 


(q; q), (=el) i4), (—( SL De 
(2 ia), C«179:9, (7 079 :4).. 


A= 


[Ref. 12: pp. 194-195] 
4. Little gq — Jacobi Polynomials 


a. Definition 
N nl 
q , abq 
| ‚9 a 
aq 


b. Orthogonality Relation 


S (bq : q). 
Y Pala, bq) r G (aq) 
x=0 GE , q)x 


(q ; q), (1—abq) (bq ; q), (abq* ; q)... (ag) à 
= —— On — (a 
(abq ; q), (1—abq°"t') (aq ; q), (aq ; q). d 
pere 0 «- gq, aq <l [Ref 12: p. 166] 
5. Big g—Jacobi Polynomials 


a. Definition 


_N n+1 
Co IA x 
sia eid ada | “| 
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b. Orthogonality Relation 


E (arg) ala ae 
P (x;a,b,c;g) P (x; a, b, e; q] ——— a— dx 
j; n "P ? Noo (OxXIC+ Doo 4 


= yg Dn (1269) (69 ; q), (aba les Dn ac)” al) , 


(abq ; q), (1—abq^"* ) (aq ; q), (eq ; q), 


m,n 


where 


- [4 0: diss Grit idis 
E | (x ; qs (bx/c , ae a 


aq (1—4) (q ; q),, (cla ; q), (aqlc ; q),, (abq? ; q),, 
(ag Q) (bq (CG; Gard ag. 


[Ref. 12: pp. 167-168] 
6. q -Krawtchouk Polynomials 


a. Definition 


i E x, Sane 
K(x;a, N;q) 930 n 0 +4, q 


b. Orthogonality Relation 


N em A 
= a , X 
) K q oe ol gm ang) m 25 (-a)” 
, x 


x=0 
(T) (q I Dn (1 T a) (ag ; q), 
2 


- (-qa ^ ;q)y a^ q ÑO TW n 
(a le CR Con 


x (-aq^ ^y" ai aie Simon 


iNet, 12) pass} 
7. q—ĦHahn Polynomials 


a. Definition 


DES abq"* q^ 
QC y A, b, N; q) = 302 | o N q q f 
aq, q 
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b. Orthogonality Relation 
| (aq ; q), (bq ; q) 
y (x;a,b, N 5 q) Q(x ; a, b, N  q) ——————À—-—- 
x=0 " | A e 


(abq’ q)w(aq) "` (q ; 4), (1—abq) (bq ; g), (abq^ ^ ; qr 
CHOR (abq ; q), (1—abq*) (aa 3). (4 ; q), 


()-n 
x (—-aq) q Oe etal =O) ses Ni 


(aq) ^ 


Ber 2: p. 165] 
8. Dual g —Hahn Polynomials 


a. Definition 


q ERE cq” N 
Rhe Nd ha Tae “| 
au Oca 


where u(x) 2 q^* teg 
b. Orthogonality Relation 


N 
Y Ral) 5 b, 0, NV 9) Rallo); b,c, N; 9) 


x=0 


(ca: q), (1—cq7*) (bea ; a). (4 ^; )x 
Ge inc ae. 
— Clergy — (G5 Mn (Odi Dn 
(q — (beqi qu (a ^ iq), 


. G) ( * 


(ca ^)" ALO 
IE 12: p. 166] 


9. g—Racah Polynomials 


a. Definition 


= ! up EUR 2 cdg**' l 
| MET, 
aq, bdq, cq 


n 
q 
(8) ni ; a, b, c, d ; q) 9 405 | 
where p(x) = q-* + cdq! 
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b. Orthogonality Relation 
When aq, bdq, or cq=q-* for a positive integer N, then 


N 
(9) » Pl u(x) ; a, b, C, d q) fo ms) , a, b, C, d; q) w(x ; a, b, C, di q) = h, OA 


x=0 
where 


w(x ; a, b, c, d ; q) 
(10) (cdq ; q)x (aq ; q), (bdg ; q). (ca ; 4)& | — cdq 


© (qq); ("cdq ; q)x (b™'cq ; q)x (dq; q)x (aba) 1— cdq 


2x41 


and 


_ (a; q), (1—abq) (bq ; q), (ad a ; a), (abc ^a ; a), (cda 
— — (abq;q), (1—abq*"*) (aq ; q), (bdq ; q), (cq ; q), 
(cag a) en ais Cine a) a en 
(a 4), (a 5 q iq), 


(11) 


(Ref. 20: p. 1014]. When ag, bag, or cq is equal to q~”, the infinite products in (11) 
reduce to finite products. Hence the orthogonality relation (9) is valid for allg provided 
no zeros are introduced into denominator terms [Ref. 18: p. 4]. 

10. Askey-Wilson Polynomials 


a. Definition 


Pax VN NET | q) =a (ab : Dn (ac Da (ad ; Dn 


(12) Gna abedg” , ae”, ae 
X 403 Ged 


ab, ac, ad 


where x= cos Ore Tepa] 
b. Orthogonality Relation 
For- l <a oee ta ale 


l w(x ; a, b, c, d | ix Š 
HE | Br las b, c,d | q) prix; a, b, C | q) Bm = h, Ò m,n 
—| 


E Js 
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where 


|] (1 — 2(2x° - 1)4* 4) 
w(x ; a, b, c, d | 1) = AGB WO AG) 
with 
OT. (1 — 2axq" 4 a^4^*) 
k= 
and 


, _ _tabeag” 5 dan (abcd ; q), (a ;4)> (aba; on 
"> (acq" ¡qdo ladg”; q), (6€9" 5 Dao (04q" 5 Doo (C49" 5 Don 


[Ref. 18: pp. 11-14] 
We can establish a formal connection between these last two classes. If the 


g —Racah polynomials are written 


à q. O ga CG e 
4223 a'q, b'd'q, c'q | q, q 


and the following parameter changes are made: 


EGG , _ ba Dod E 
a' = RE r Z 


(or likewise using any permutation of { 6,c,d} assuming at least one of these is non- 


zero), as well as the change of variable q~ = ae”, then we obtain 


a”, abedq" |, ae”, ae” 
403 q q 


ab, ac, ad 


with 


u(x) um UE T CONGU 


7 E 
= ae” +ae TE 2a cos 0 


9n 


That is, the 4 — Racah and Askey-Wilson polynomials are virtually the same, differing 
only in their parameters, normalization, and variable. 

Because the g—Racah and Askey-Wilson polynomuals are essentially the 
same, the names are often used interchangeably in the literature. We have presented 
these polynomials as distinct classes in order to emphasize the continuous and discrete 
natures as expressed in the orthogonality relations. 

As stated at the beginning of Chapter IV, the “classical orthogonal 
polynomials” are defined to be those which are special or limiting cases of the Askey- 
Wilson (21) or g—Racah polynomials (18) [Ref. 31: p. 57. We can now make this 


statement a bit more precise. Letting 


in (9) and (10) and taking the limit as q—> l, the ordinary Racah polynomials 
R( 2(x):a, B, y, 6) and their weight function as defined in Formulas (30) and (32), 
Chapter IV, Section C.3 are retrieved. A similar limiting process will recover any ordis 
nary orthogonal polynomial class from its q —extension. Moreover, any q —extension 
previously discussed is a special or limiting case of the q —Racah. For example, letting 
cq = q? and d 2 0 produces the q —Hahn, etc. In this way, we see that the self-dual 
g —Racah polynomials encompass all the previous classes. They also satisfy three-term 
recurrence relations, second order difference equations, and Rodrigues’ formulas with 
respect to g —divided difference operators. The interested reader is referred to [Refs. 18, 


12, 20] for details about these very rich classes. 


C. CONCLUDING REMARKS 

In Chapter V, we presented a few of the traditional applications in which orthogonal 
polynomials have played an important part. More recently, an enriching mterplay has 
developed between the theory of orthogonal polynonuals and other mathematical and 
mathematically-related areas. 

Efficient computational methods have been devised for determining the many useful 
quantities associated with orthogonal polynomials (such as their zeros, recurrence coef- 
ficients, etc.) [Ref. 32: pp. 181-216]. Various classes of orthogonal polynomials have 
also played a role in digital signal processing [Ref. 32: pp. 115-133], quantum mechanics 
[Ref. 32: pp. 217-228], and birth, death processes (Ref. 32: pp. 229-255), (Ads ances mere 


field of combinatorics and gruph theory have allowed new geometric interpretations of 
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orthogonal polynomial identities, some of which have very important consequences for 
“association schemes” and the designs of codes [Refs. 33, 16, 34, 32 : pp. 25-53, 35]. 

Physicists have introduced various versions of “diagrammatic methods”: ways of 
understanding orthogonal polynomials through 3-j and 6-j symbols, their generaliza- 
tions, and accompanying identities by formally associating them with pictorial schemat- 
ics representing forces of physical systems that conserve angular momentum [Ref. 36]. 
Powerful new techniques involving “quantum groups" have been used to generate new 
identities for some classes [Ref. 32: pp. 257-292]. Further investigation into the 
electrostatics problem discussed in the text has led to the formation of the famed 
"Selberg beta integral” and its generalizations. It has yielded to analysis via the study 
of the “root systems” of Lie algebras, and has found applications ranging from statistical 
mechanics to computer algorithm complexity [Refs. 30: pp. 48-52, 32 : pp. 311-318, 37]. 

Finally, research into the general structure of q —series has led to many surprising 
connections, and is intimately related to the many remarkable and powerful number- 
theoretic formulas discovered by S. Ramanujan, the famed Indian mathematical prodigy 
[REIS SO: pp. 87-93, 38: pp. 55-66, 39]. 

This utility and promise of future applications provide ample motivation and justi- 


fication for continued study of the intrinsic structure of orthogonal polynomuals. 
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